A direct transition between a Neel ordered Mott insulator and a d x 2_ y 2 
superconductor on the square lattice 
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In this paper we study a bandwidth-controlled direct, continuous, phase transition from a Mott 
insulator, with easy plane Neel order, to a fully gapped d x 2_ y 2 superconductor with a doubled unit 
cell on the square lattice, a transition that is forbidden according to the Landau paradigm. This 
transition is made possible because the vortices of the antiferromagnet are charged and the vortices 
of the superconductor carry spins. These nontrivial vortex quantum numbers arise because the 
ordered phases are intimately related to a topological band insulator. We describe the lattice model 
as well as the effective field theory. 
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One of the central questions in the theory of the high 
T c cuprates is how the antiferromagnetic (AF) insula- 
tor is converted into a d-wave superconductor (SC) upon 
doping. If a continuous transition between these, in some 
sense, diametrically opposite phases could be realized, it 
would be possible to construct an unified description of 
the SC and AF in terms of a single critical theory^. 

An immediate indication that such a theory would 
be rather unconventional, is that it violates the 'Lan- 
dau' rules which dictates that continuous phase tran- 
sitions are only possible when the symmetry groups of 
the two phases have a subgroup relationship. It has 
been recently proposed that Landau forbidden transi- 
tions could in principle occur in quantum systems^, via 
'deconfined' quantum critical points. In this area, early 
work on the AF to valence bond solid transition on the 
square lattice 2 -! 3 , has been generalized to a number of 
other transitions^j&i. However, except an early numer- 
ical claim£, the physically most interesting case of a di- 
rect transition between an AF insulator and a SC has 
not been discussed before. Here, we provide a theory for 
such a transition. 

The essential physics here is that vortices of the AF 
(in our case the easy-plane variety) are charged and the 
vortices of the SC carry spin. The condensation of either 
type of vortices drives the system between the two phases. 
The reason that the topological defects of these two 
rather conventional phases carry unconventional quan- 
tum numbers is because both phases are intricately re- 
lated to a topological band insulator. The later is a spe- 
cial type of band insulator which can be differentiated 
from conventional ones by a topological quantum num- 
ber. 

The field theory that emerges at criticality is the 
easy plane version of the non-compact CPi (NCCPi) 
model 2 ^. An important theoretical feature is that in con- 
trast to most earlier theories 4 -^ there is no emergent U(l) 
symmetry at criticality. An interesting critical point be- 
tween an magnetic insulator and a d-wave SC with nodal 
quasiparticles was proposed in Ref.— . However, in con- 
trast to the present theory, it is hard to establish that 



the insulating state is truly an antiferromagnet. 

In contrast to the cuprates, where the transition is trig- 
gered by doping, here we consider a bandwidth controlled 
Mott transition, where the system remains at half filling 
throughout. Moreover, in addition to the usual d-wave 
symmetry, our SC breaks lattice translation and time re- 
versal symmetry due to a staggered spin dependent hop- 
ping term. This leads to a gap appearing at the nodes. 
Doping triggers a gap closing transition and gap nodes 
are restored at sufficiently high doping. Finally, our AF 
phase has a spin anisotropy favoring easy plane magnetic 
order. Topological defects with nontrivial quantum num- 
bers were also discussed in ReL 10 ! 11 ' 12 , in the context of 
cuprate superconductors. 

While in mean-field theory the easy- plane NCCPi 
model clearly predicts a continuous transition, the effect 
of fluctuations has been actively debated ^ 13 ! 14 ' 15 ! 16 ' 17 . 
If we consider the general family NCCP„ models, these 
are know to have continuous transitions for n = and 
n — > oo. For the case of n = 1, with full SU(2) symme- 
try, the case for a continuous transition has been made in 
some microscopic realization a 14 i 15 ' 16 , although other au- 
thors have concluded that weak first order transitions are 
present in the entire parameter ranged 7 -. Since the nature 
of the transition depends on the microscopic realization, 
it is clear that more numerical studies on a variety of 
models, including those with easy plane symmetry, are 
required to settle this question conclusively. In this pa- 
per, we do not further address this issue, rather we aim 
at establishing whether a continuous phase transition be- 
tween the d-wave SC and AF insulator that is possible 
even in principle. 

The plan of the paper is as follows. We start with 
the Bogoliubov-de Gennes description of a d x 2_ y 2 super- 
conductor with a particular unit cell doubling shown in 
Fig. [TJ (see later) which leads to gapped nodes. It turns 
out that this is a topological superconductor, which in 
the presence of particle-hole symmetry, has stable edge 
modes. We show that vortices of this superconductor 
carry spin and are bosons. We ask what happen if these 
vortices condense. This is answered by studying the field 
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FIG. 1: (a)The arrows indicate the direction in which the 
nearest neighbor hopping amplitude is e 1 ^ (hopping in the 
reverse direction has amplitude e~ 1 ^. The color indicate the 
sign of the next nearest neighbor hopping, Red means nega- 
tive. (b)Bulk and edge spectra of Eq. ([T]) and Eq. (2). We 
used parameters <f> — n/5 and t = 0.2 (see the discussions 
below Eq. (JTJ). 



theory of the vortex excitations, which predicts that the 
vortex condensed phase is an easy-plane antiferromagnet 
and the critical theory is an easy plane NCCPi model. 
A duality transformation allows us to study the transi- 
tion in reverse, where condensation of charged vortices 
of the AF order leads to superconductivity. The self du- 
ality of the NCCPi modelS. leads to an interesting spin- 
charge duality, whose consequences are described in Sec- 
tion IIV Al As a consistency check, we study the same 
problem from the insulating side. Quite interestingly the 
most natural slave boson formulation leads to the same 
phase diagram and critical theory. This alternative ap- 
proach also clarifies the role of symmetries required to 
realize this transition, and allows us to generalize this to 
other lattices. In the conclusion we discuss the similarity 
of the phenomenon described in this work and that of the 
cuprate superconductors. Eight appendices contain the 
important technical details of various part of the paper, 
including Appendix [F] where the nature of the magnetic 
order are studied in details and Appendix [H] describes 
an equivalent theory on the honeycomb lattice, where a 
transition between an easy plane AF and a triplet super- 
conductor occurs. 



I. THE TOPOLOGICAL D-WAVE 
SUPERCONDUCTOR 

Let us begin with a mean-field description of one side 
of the phase transition: the d-wave superconductor. The 
usual d-wave superconductor has gapless nodal quasi- 
particles; here we introduce a modification which leads 
to a unit cell doubling, and a full gap around the Fermi 
surface. The mean-field Hamiltonian describing this su- 
perconductor is given by 



H sc _ ^2 Re(xij)c} (T Cj (7 - E r) i Im(xij)e<r<r'Ci <T c 

(ij)<7 (ij)crcr' 

+ E crfijcLcjV + h.c. 



(1) 



Here i,j run through the sites of a square lattice, r\i is 
the AF staggered sign, and denote the near- 

est and second neighbors, respectively; Xij = e±l ^ an d 
Uj = ±t (t — real). We note, the nearest neighbor hop- 
ping, Re(xij), is the same as the hopping parameter asso- 
ciated with the "staggered flux phase" of the Heisenberg 
antiferromagneljiS, (for a review, see RefJ^). We illustrate 
the Xij an d Uj patterns in Fig(QJa)). It is important to 
note that the pairing term has d x 2_ y 2 symmetry; due 
to the presence of r\i in the second term of Eq. (JlJ, the 
pairing is translation-invariant. 

However, because of the spin-dependent second neigh- 
bor hopping, both lattice translation T and time reversal 
TTZ are broken. However, T o TTZ and reflection along 
the vertical(horizontal) line passing the center of each 
plaquette (V) are symmetries of the Hamiltonian. ((By 
looking at Fig(QJa)) one might think that reflection is 
not a symmetry. However this figure merely illustrates 
the pattern of Xij and Uj not the symmetry of the Hamil- 
tonian.) In addition, Eq. |T]) has spin S z rotation sym- 
metry, and the charge U(l) symmetry is, of course, spon- 
taneously broken. 

We will see in later discussions that T o TTZ and 
V together with the U(l) — S z conservation and the 
[/(l)-charge conservation (which is spontaneously bro- 
ken in the current superconductor phase), ensures that 
the direct SC<->AF transition can be continuous at half- 
filling. There are two more symmetries of the Hamil- 
tonian Eq. ([I]) which simplify our analysis but are not 
necessary for the direct transition to occur: the particle 
hole conjugation PH (c^ — > Jjic], , Cj| — * —T]iC^) and and 
translation together with Ci a <-> Ci_ CT (To "f<->J,). 

Eq. |T]) is closely related to a special type of topological 
insulator (TBI) Hamiltonian 

""" E xvfUs*+ E rtijfifj* +h.c..(2) 

Indeed, one can transform Eq. |T]) into Eq. J2]| by making 
the following transformation 



V2 



Cil 



f, 



a 



V2 



(3) 



If we drop the a in front of the Uji Eq. ^ turns into 
the Hamiltonian for the "chiral spin liquid"—. With the 
<7, the chirality of the spin up and spin down fermions 
are opposite. In that case Eq. §2$ describes a band in- 
sulator with opposite "TKNN index"— (±1) for the spin 
up/down fermions. Hence this type of topological insu- 
lator will exhibit the spin hall effect. However, unlike 
the spin hall insulator introduced in Ref.— , our spin hall 
insulator breaks time reversal symmetry. (To appreciate 
the full relationship between Eq. (fTJ) and Eq. @, see Sec. 

m 

Since Eq. (|3|) is an unitary transformation the excita- 
tion spectra of Eq. ([2|) and Eq. (p} are the same. One 
only needs to reinterpret the electron/hole excitations in 
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FIG. 2: The Bogoliubov quasiparticle (q.p.) spectra of a pair 
vortex and anti vortex; The red lines mark the core states, 
and the black solid line represents the Fermi energy. In the 
limit of infinite vortex-antivortex separation each core level is 
doubly degenerate. We used parameters <j> — n/5 and t — 0.2 
(see the discussions below Eq. (HJ). 

the TBI as the Bogoliubov quasiparticle excitations in 
the superconductor. Interestingly, the conserved S z of 
Eq. @ translates into conserved S z of Eq. fl} as well. 
Among other things the above correspondence implies 
that the superconductor described by Eq. (Q]) has gap- 
less edge excitations (FigfTJb)). Thus it is justified to 
regard this type of superconductor as "topological". In 
Appendix IG 21 and near the end of Sec lVBI we show that 
the edge modes requires the PH symmetry which is diffi- 
cult to realize in a material. Therefore in general there is 
no gapless edge modes and the SC phase is not topologi- 
cal. On the other hand the physics of the direct SC^AF 
transition, which is the main result of the current paper, 
does nog require the existence of gapless edge modes. 

Now let us study the vortices of the superconductor. 
Direct numerical calculations for a vortex and antivor- 
tex pair gives the spectra illustrated in Fig. (|2]) , where 
each core level is doubly degenerate (in the limit of infi- 
nite separation between the vortices) / There AV and BV 
stand for a vortex located on A dual sublattice and an 
antivortex on B dual sublattice, and etc. Here A/B dual 
sublattice locate at the center of plaquette in Fig.(QJa)) 
marked by red/black cross lines. The core levels of these 
two types of vortices lie on opposite side of zero. In the 
ground state all levels below/above zero are filled/empty. 
For each vortex/antivortex, there are three low lying ex- 
cited states arise from different way of occupying the two 
core levels associated with it. Explicit calculation of the 
ground state value of S z associated with the four panels 
of Fig.@ gives 

SXv = Slv = 1/2, and S^ v = S*# = -1/2. (4) 

The spin distribution associated with (AV,AV) and 
(AV,BV) is shown in Fig.©. 

The fact that vortices carry spins is the most interesting 
aspect of the superconductor under consideration. 

In the above numerical study, the size of the vortex core 
(the region in which the pairing amplitude is suppressed 
to zero) is one lattice constant. As the core size of the 



FIG. 3: The S z distribution of an (AV,AV) and (AV,BV) 
pair. The total integrated values of S z around each vortex is 
5 Z = ±1/2. 

vortices increase, their core levels move toward zero. As 
the core size diverge, a continuum BdG equation applies. 
In that limit all core levels are situated at zero energy 
(Appendix |B|) and the previous three (low lying) excited 
states become degenerate with the ground state. When 
that happen it is no longer possible to differentiate the 
spectrum of A and B dual sublattice vortices. In this 
limit each vortex/antivortex has four degenerate states: 
two with S z — ±1/2 and the other two with spin zero. 

Next we determine the statistics of the vortices. For 
this purpose it is important to note that under 

PH: AVf <-► AV T ,BV X <-> BVj 
T o T1Z: AVf «-> BV X , BVj <-> AV T (5) 

Now let ^BVi,BVi ! and etc., be the phase factor due to the 
(counterclockwise) exchange of two BVj, vortices. ToTIZ 
plus PH imply </>BVj.,BV; = <^av av t • On ^ ne other hand, 
To f«-»J, requires 0BV i ,BV i = 0av t ,av t ■ Put together the 
above implies </>av t av t = 0av t ,av t - Consequently 

0AV T ,AV T = ±1, (6) 

i.e., AVf is either boson or fermion (the same for BVj). 
Although we used the symmetry To "f<->J, to derive this 
result, the result Eq. J6]) does not depend on that symme- 
try because statistics cannot change so long as the core 
levels do not merge into the continuum. 

To further distinguish the two we consider the follow- 
ing gauge transformation a a — ► e 10 *' 2 ^. This transfor- 
mation remove the vorticity in the paring order param- 
eter up to a 7r cut. The hopping parameters acquires 
a 7r cut as well as a smooth hopping phase e laij where 
otij = Mod( 6 ' 2 6j , 7r). This hopping phase is often re- 
ferred to as the "Doppler shift term" in the theory of su- 
perconductivity. If one gets rid of the hopping phase, and 
leaves the 7r cut in both the pairing and hopping param- 
eters, both core levels of a given vortex/antivortex move 
to zero energy and there are four degenerate vortex states 
, like the case of the diverging core size discussed earlier. 
In Fig. ([4]) we illustrate the evolution of the spectra for a 
AVf and a BVj vortex after the following modification 
of the hopping phase: 

e ia « -> e iXa ^ < A < 1. (7) 
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FIG. 4: The spectra of a AVf and a BV| vortex as a function 
of A. From left to right A = 1, | and 0. 



In the A = limit Eq. transforms Eq. dXJ) to a spin- 
Hall insulator with Zi flux in the plaquette previously 
occupied by the vortices. This is a problem we have 
previously solved and the S z = ±1/2 vortices (which 
are referred to as spin fluxons in Ref^) are bosons^. 
It is simple to show that Eq. © holds true for all A 
hence Eq. ^ remains valid. This, plus the fact that the 
spin fluxons adiabatically evolve into the spin-carrying 
vortices, implies 



W T ,AV T 



1. 



A=l 



(8) 



, i.e., the AV| vortices are bosons (same argument applies 
to BV|, AV t and BV |.) 



II. THE SUPERCONDUCTOR- VORTEX FIELD 
THEORY 



Since the vortices are bosons, they can condense. The 
condensation will trigger the destruction of superconduc- 
tivity and lead to an insulator. However, since the vor- 
tices carry spin, we expect the resulting insulator to have 
magnetic signature. To understand what happens after 
vortices have condensed, we examine the effective field 
theory describing the phase fluctuations of the super- 
conductor. Phase fluctuations appears in two different 
types: gaussian (spin wave ) type and vortex-antivortex 
fluctuations. The latter is responsible for destroying the 
superconducting state when it is strong. For example, it 
is the vortex condensation that triggers the ordinary su- 
perconductor to insulator transition. A convenient way 
of describing the vortex physics is to use the dual form 
of the phase fluctuation theor y 24 i 25 . In this form, the 
vortices are particles and the original boson density and 
current become magnetic and electric fields respectively. 
Since this dual form is widely used in literature, we will 
refer the readers to reference s 24 ' 25 for details. The rele- 
vant field theory for the superconducting vortices in three 
dimensional Euclidean space takes the form: 



■«I*tI 2 I*J 



f (V x a) 2 + fi-(Vx a\9) 



Where 'fj,^ are relativistic vortex fields, and for 
brevity, we have omitted the space-time dependence of 
all fields. In addition, we have used a three vector nota- 
tion V = (do, A^, A y ) to denote space time derivatives. 
The connection between these relativistic fields and the 
microscopic vortex fields is discussed in Appendix [A} As 
usual, a gauge field a is introduced such that x a — 
is the Cooper pair 3-current caused by the gaussian fluc- 
tuations of the pair-condensate phase. Hence, it is cou- 
pled to the external electromagnetic field A, contained 
in the last term of Eq. The minimum coupling 

Y*n \( 9 m + ia v - i §^°' 3 )*<r| 2 between and re- 
flects the fact that vortices see the 3-current of the con- 
densate as electric and magnetic fields. To probe the spin 
response we have we coupled a "spin gauge field" to 
the conserved S z current. Note that both A^ and A^ are 
non-dynamical. The fact that the vortices minimally cou- 
ple to the spin gauge field reflects the fact that they carry 
spin. The |tr 3 in the coupling between ^ and A^ re- 
flects the fact that the boson carries S z = ±1/2. Fi- 
nally, we have included standard Ginzburg-Landau type 
quartic ^ a terms. They describe the interaction between 



the order parameters of the two type of vortices. 



To ensure the correctness of the low energy effective 
theory Eq.®, the following conditions must be met: (1) 
vj/| — > vj^ symmetry (the two boson fields are inter- 
changeable), (2) $ — > vfrf symmetry (particle and anti- 
particle are on the same footing which requires the ef- 
fective theory to be relativistic), (3) ^ particles do not 
experience a background magnetic field; i.e. there is no 



linear term h(d x a y 



-d y a x 



in the action above. We show 



in the following that the three requirements are ensured 
by the symmetry operations PH, V and ToTlZ. (Even- 
tually, in Sec lVIIl we will show that even PH is not re- 
quired: as long as the density is fixed at half-filling, the 
low energy effective theory Eq.@ is correct.) 



Note that according to Eq.© both PH and V trans- 
forms a vortex into an anti-vortex, while preserving spin 
while ToTlZ transforms a vortex into an anti- vortex and 
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flips spin. Therefore, 

PH: => ^t-f^.^-f*} (10) 

V- => *iT->4(01'*l->4 W t ( n ) 
ToTIZ: => * T -> -> (anti-unitary) 

(12) 

Although P and Pi/ perform similar transformations 
on the vortex field, they act very differently on the back- 
ground flux (d x a y — dya x ). V changes the sign of this 
term while PH leaves it invariant. Therefore the com- 
bination of V and PH ensures condition (3). Similarly, 
Eq. {TP} and Eq. (fT2"j) validate the conditions (1) and 
(2). Therefore in the presence of Pi/, P and T o T1Z 
conditions (1),(2) and (3) are satisfied. This dictates the 
effective field theory above. In later discussions we will 
derive the same theory from a completely different ap- 
proach Eq. ([39]) . This will clarify the role of particle hole 
symmetry, which will not be necessary as a microscopic 
symmetry to achieve this transition. Note, in the absence 
of the external fields A, A s and in the presence of easy 
plane anisotropy (u < 0), the above theory is identical to 
the easy plane NCCPi of Ref£. 



III. THE AF ORDER IN THE VORTEX 
CONDENSED PHASE 

Let us now focus on Eq. ©; for m 2 large and positive, 
the 'J'o- vortices are absent at long wavelength/ time. In 
that limit Eq. becomes a gaussian theory in a M . Inte- 
grating out a M generates a Higgs effective action for A^, 
in other words, the system is a superconductor. When 
m 2 changes sign the spin-carrying vortices will condense. 
If u > the and ^ vortices suppress each other, 
hence energetically it is favorable for (^j ) ^ 0, (^i) = 
or (V^f) = 0, (t&i) ^ in the condensed phase. When 
-4« < u < it is favorable for (1) (* T ) = (vf^) ^ 
or (2) ^ while (* T ) = (%) = 0. In case (2) 

the order parameter has no vorticity (hence is an local 
operator) while flips the spin. When it develops expec- 
tation value, superconductivity needs not be destroyed. 
Since (*{*;) ~ (S + ) ^ 0, the system exhibits SC 
and magnetic order simultaneously. While this scenario 
can certainly be realized in certain parameter regime, in 
the rest of the paper we focus on the more interesting 
case where the system goes from = ("J^) = to 

= (Wj) 7^ in a single transition. 

When single vortices condense, a SC turns into an in- 
sulator. When both and are none zero, (^j^) 
hence (S + ) are also non vanishing. As a result, the in- 
sulating phase has easy plane magnetic order. However, 
what is the order pattern? Is it AF or ferromagnetic? 

To minimize technical complexity, here we present a 
simple, but less rigorous, way to clarify the nature of the 
magnetic order. We put the more rigorous calculation, 
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FIG. 5: For Re(x) = 1, t = 0.8, and Im(x) = 0.3 (see Eq.©) 
and a vortex located in the plaquette labeled by ©, we plot 
the magnetization (Si) on site i in a small in-plane field B — 
O.ln. The sign of the magnetization is labeled by the arrows, 
and the length of the arrows is proportional to The 
calculation is performed on a 1152-site system with one vortex 
and one anti- vortex separated by half system size. 

based on computing vortex tunneling amplitude, in Ap- 
pendixlEl To mimic the coherent tunneling (('J'j'i'j) ^ 0) 
between the spin up and spin down vortex states in the 
magnetic phase, we turn on a small in plane magnetic 
field (the field strength needs to be larger than the spac- 
ing between the core levels in the first panel of Fig ^ to 
mix the spin-up and spin-down vortex states in a single 
static vortex. In FigfS] we show the expectation value of 
the in-plane spin component in the vicinity of the vor- 
tex core when there is an in plane magnetic field in h 
direction. The AF order is clearly seen. 



IV. TRANSITION FROM THE AF SIDE 

In the previous section we have analyzed the SC <-> 
AF transition from SC side using the effective action 
given in Eq.®. In this section we approach the tran- 
sition from the AF side, where the AF order is destroyed 
by the condensation of AF vortices. Since superconduc- 
tivity emerges after this condensation, it suggests that 
the AF vortices carry electric charge. In the following we 
show that this is indeed the case. 

Let us assume in vortex condensed phase the U(l) or- 
der parameter (vPf) and (^j.) have equal magnitude. In 
the following we derive an effective theory describing the 
vortices of these order parameter fields. The AF— +SC 
transition is then triggered by the condensation of these 
new vortices. The derivation is accomplished by the per- 
forming duality transformatio n 24 ! 25 on Eq. © (for details 
see appendix fA")). In fact, in the absence of the external 
gauge potentials A, A s , this has been carried out in Rcf. 9 
where it was discovered that the theory at criticality is 
self-dual. In appendix [A] we extend it to the case with 
these potentials present. The end result is the effective 
theory given in Eq. (|13[) . Essentially in the AF phase 
described by (* T ) = y/pe^ , * x = y/pe^i , the re are 
two types of vortices. One causes the winding of <fi u a by 
2ir and the other the winding of <f>\ by — 2-7T (here is the 
minus sign is chosen for convenience) . In the presence of 
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such vortices, a nucleates ±7r fluxes to screen the vortic- 
ity (so that the kinetic energy is minimized). From the 
term i^A ■ V x a in Eq. ((9]) we find that these two types 
of screened vortices carry charge ±e. We label these vor- 
tices by z\ and z 2 respectively. This explains the oppo- 
site sign in the minimal coupling to A in the first term 
of Eq. (fT3|) . The curl of the a M , namely V x a, describes 
the spin 3-current caused by the gaussian fluctuation in 



the phase of (S~) ~ e ^t-0i). This is reflected in the 
last term of Eq. ([13")) . Since both zi and z 2 creates 2ir 
winding in 0f — <f>^ it couples to a M minimally as shown 

by the first term of Eq. (JT3"| . Finally z\z2 inserts charge 
2e while has no spin vorticity (hence is an local operator) 
is the Cooper pair operator. 



dual 



E 



a=l 

ih 

2^ 



E \(d^ + ia^ - ieA^a 3 )z a \ 2 + m 2 \z a \ 2 



(V x a) ■ A" 



u\z 1 \ 2 \z 2 \ 2 + 



(V x a) 2 +ci(V x A) 2 
(13) 



Note, the form of Eq. (fT3|) is almost identical to Eq. © 
except for the exchange <-* z a and 2eA <-* HA S . Since 
Eq. (fT3")) is the dual theory of Eq.©, it must be consis- 
tent with microscopic symmetries of Eq. (J9j) . In Appendix 
IG II we further directly checked the consistency between 
microscopic symmetries and Eq. (|13l) . 

Next let us how to describe the magnetically ordered 
phase using Eq. (fl3]) . In the ^ phase, m 2 > 

and the z a bosons are absent at long wavelength/time. 
In that limit we again have a gaussian theory 

£gau g c-> 5(Vxo) 2 -^A^Vxo. (14) 
Z Ztt 

describing a "spin superfluid" . The above action is dual 
to the following Goldstone action 

£ gs = ^E(^X-^) 2 - (15) 

Here the phase angle x transforms under a angle 
(3-S z rotation as x ~ * X + P- It can also be 
shown that computing the spin-spin correlation function 
( e *x(xi:t) e - l x(x2,t)^ using action Eq. (fT5|) is equivalent to 
computing the monopole-anti monopole correlation func- 
tion using Eq. ([14"]) , i.e., 

(e ^(x llt ) e -ix(x 2l *) )£Gs = (y+( Xl , <)F -( X2;t))w .(i6) 

Here V + and V~ are the creation operators of the 
monopole and antimonopolc in a-gauge field. Since a 
monopole (anti-monopole) is the event at which ^- (V x 
a)o = S z /h changes by 1 (—1), we conclude that 

(V+( Xl ,t)V-(^ 2 ,t)) ~ (5+(xi,t)S-(x 2 ,t)). (17) 

Thus the S z spin superfluid described by Eq. (jT^J) and 
Eq. P^|) is an XY ordered magnet (from the previous 
discussion we know this order is AF) and z\,z\ create 
charged XY vortices (anti- vortices). Finally z\z 2 inserts 



charge 2e hence creates a monopole in a-gauge field, and 
^j^E'l raise S z by 1 hence is a monopole operator 
in a-gauge field. (We will come back to this point later 
in the slave-rotor description SecfVland Scc[Fl where we 
explicitly compute the monopole quantum number for V' 
and find it is AF.) 

Since Z\ and z 2 carry the vorticity of magnetic order 
parameter their condensation will destroy the AF order. 
However since z\ and z 2 carry charge, the condensate 
is a superconductor. Now we have gone through a full 
circle. We start from the d-wave superconductor, the 
spin-carrying vortex condensation brings the system to 
a XY ordered antiferromagnet. In the reverse direction, 
start from the antiferromagnet, the condensation of the 
charge-carrying vortices brings the system back to the 
superconductor. 

In Lapp lDlEl we perform the parallel study of the 
AF vortex core states as what we did for SC vortex in 
Sec lIIIIIinl We indeed find that there are two mid-gap 
levels in the AF vortex core. The filling/unfilling of these 
levels lead to charge-±l AF vortices, whose condensation 
triggers the SC order. Analogous to the discussion in sec- 
tion III we can determine the SC pairing symmetry by 
studying the tunneling property of the magnetic vortices. 



A. Spin-Charge Duality 

As mentioned previously, the self duality of the easy 
plane NCCPi model at criticality, leads to spin-charge 
duality here. The effective theories near the transition 
point, Eq. © and Eq. (13]), map onto one another with 
the substitution ^ a <-> z a ,2eA <-* HA S which implies 
that the roles of spin and charge are identical at critical- 
ity. We point out some consequences of this duality in 
the following. 

At criticality we expect both the electric and spin con- 
ductivity to be finite and nonzero. The self duality en- 
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sures the electric conductivity (cr c ) and spin conductivity 
(a a ) to be related by 2f = • g enera lj se lf duality 

allows another term in the response action : -A-Vx A s , 



i.e., a spin Hall term. However this is forbidden by parti- 
cle hole symmetry PH: A — ► —A.) Therefore we expect 
the critical gauge action to have the following form: 



-eff 



2 



(V x eA),. 



= (V x eA\ 



(V x -A s ),- 



7 V ; 



h 



= (V x ^A S ) M 



(18) 



Note, in a regular Fermi liquid, a Wiedemann-Franz like 
relation connects the metallic spin and charge conductiv- 
ity in precisely the same way, since at low temperatures 
both quantities are transported by electronic quasipar- 
ticles with a fixed charge to spin ratio. At the critical 
point described above, although we have spin-charge sep- 
arated excitations, the self duality of the theory restores 
this Wiedemann-Franz like relation between the metallic 
spin and charge conductivities. 

At criticality the correlations of the magnetic order 

I/2+77/2 

parameter < S + (r,t)S (0,0) >~ r *_Ji t i and 
the Cooper pair order parameter < A*(r, t)A(0, 0) >^ 
_^_^_^i/2+i?/2 related by the self duality, and con- 
strained to fall off with the same power law. Since these 
are both expressed as bilinears of vortex operators, the 77 
exponent is expected to be large. Thus, the critical point 
unites the antiferromagnetism and superconductivity in 
a remarkable symmetric fashion. 



V. SLAVE BOSON FORMULATION 



In SecU we started from a lattice Hamiltonian for 
a superconductor and showed that the vortices carry 
spin. We then condensed these vortices and showed 
that, as a result, the system becomes an AF ordered 
Mott insulator. In this section we begin with a lattice 
theory for this Mott insulator, and show that we can 
arrive at the same conclusions from this different angle. 
The slave boson gauge theory approach used below, cap- 
tures both the superconducting and antiferromagnetic 
insulator phases within a single formulation. While the 
former is obtained by the Higgs mechanism, the latter, 
surprisingly, is the Coulomb phase of the gauge theory. 
This is an interesting counterexample to the standard 
folklore that slave boson theory can readily describe 
superconducting states but not magnetically ordered 
phases. Here we establish that magnetic order of the 
easy plane variety can indeed be captured. 



A. Insulating State 

Consider a half filled strongly correlated electronic 
model on square lattice. 

#Hubbard = ffhop + t^K-l) 2 (19) 

i 

where the first term is a hopping Hamiltonian, while the 
interaction term tends to prefer single occupancy per site 
(rii is the electron number on site i). Unlike most other 
such models which exhibit spin SU(2) symmetry we as- 
sume that there is only S z spin rotation symmetry;this 
is incorporated in -ffhop via e.g. spin dependent hopping 
terms. 

Deep in the Mott insulating state, when we turn up 
the value of U, the physics is described by a spin 1/2 
Hamiltonian. For concreteness we consider obtaining the 
following simple Hamiltonian: 

H = J S r ■ S r > + flatus (20) 

(rr 1 ) 

flanis — ^~ ^ ^ (S^ S r /* -\- S r Sy//), (21) 
««•"» 

where the SU(2) breaking term comes in as a ferromag- 
netic second neighbor interaction between planar compo- 
nents of the spin. The ground state of this Hamiltonian 
is expected to be an easy-plane Neel ordered state with 
spins in the XY plane. 

Below we will show that the same conclusion can be 
reached via the slave boson formulation in a controllable 
fashion. This new mechanism for obtaining magnetically 
ordered states in the slave boson approach is not just of 
academic interest. It allows us to simultaneously cap- 
ture the magnetically ordered states and superconduc- 
tivity (the latter is readily realized using slave bosons as 
the Higgs particle) within the same formulation. This 
is achieved without having to resort to 'confinement' 
which in the slave boson context does not allow for a 
unique identification of the ordered state. In contrast, 
the approach presented here will allow us to identify both 
the XY antiferromagnet and the superconducting phases 
uniquely. Moreover unifying these two phases within a 
single formulation leads to the possibility of a direct con- 
tinuous transition between these two phases. The slave 
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boson approach also allow us to identify the nature of 
the superconducting state and the symmetries required 
to accomplish such a transition. 

Consider treating the Hamiltonian (j2"0|) using a 
'Schwinger fermion' representation of spins, i.e. intro- 
duce two component fermions (/rT, frl) a t each site with 
the constraint: 

E/rUW = l (22) 

imposed on every site. Then, the spin 1/2 opera- 
tors are S r — \fra&a,a' fro' ■ Substituting this in (|2"01) 
leads to a Hamiltonian quartic in the fermionic op- 
erators. To make progress, consider the mean field 
approximatio n 18 i 26 i 27 , obtained by replacing (/J CT / r 'o-') = 
Xrr'.a.a' and treating the constraint on average. It is well 
known that in the absence of the spin anisotropy term, 
the optimal mean field ansatz with uniform amplitude is 
the staggered flux state 

H MF = - E Xrr' fUr'v (23) 
(rr') 

with the nearest neighbor hopping amplitudes \rr' = 
|^| e ±»4>/4 a i on g (against) the directions shown in the Fig- 
ure [lla) • The optimal variational state within this set is 
reached for <E> w 0.347Ti2S 

Now let us try to minimize the mean-field energy due 
to the spin anisotropy term 

=y E (f^fr"M"lM)- (24) 

«"-"» 

One way in which this can be achieved is if the two ex- 
pectation values are real but with opposite signs, i.e. 

H MF = - E Xrr'fLfr'a - E at rr" fra fr" a (25) 

(rr') <(rr"» 

the spatial structure of the spin dependent next- neighbor 
hopping t" ,i is shown in the figure [TJa). This choice of 
t s rr ,i opens a gap at the Dirac points and hence leads to 
a better mean field energy. 

The mean field Hamiltonian above has a familiar form. 
In the absence of spin dependence of the next neighbor 
hopping, this is the "chiral spin liquid" Hamiltonian 20 . 
The spin dependence leads to opposite Chern numbers 2 ^ 
(±1) for the spin up and spin down fermions, and 
hence describes a spin Hall insulator. However, unlike 
the honeycomb lattice spin hall insulator introduced 
in Refi22, this band structure breaks time reversal 
symmetry, although time reversal invariance combined 
with a translation is still a symmetry. Note also, that the 
unit cell is doubled once these spin dependent hoppings 
are introduced (the staggered flux mean field theory 
in contrast, is effectively a translationally symmetric 
ansatz). 



Gauge fluctuations and Neel order 

The U(l) symmetry, corresponding to the conservation 
of the / fermion number, of Eq. (|2"5|) is not a physical 
symmetry, but rather corresponds to a U(l) gauge re- 
dundancy. Including the corresponding gauge field a rr > 
on bonds of the lattice, requires \rr' — > Xrr'e mrr ' and 
t s rr „ — > t s rr „e larr " . The temporal component of the gauge 
field clq arises as the Lagrange multiplier imposing the 
constraint (|2"2"|) . 

Since the fermions are gapped, they can be safely in- 
tegrated out to obtain an effective action for the gauge 
degrees of freedom. In order to interpret this action it 
is useful to introduce an external 'spin' gauge field Af r , . 
This is possible since S z is a conserved spin density. Re- 
placing a rr > — > a rr > + T}Af r ,a z : and integrating out the 
fermions, one obtains the Euclidean Lagrangian: 

C E [a,A s ] =jA/.Vxa+5(Vxaf + (Vx ^A s ) 2 } 

where we have used a three vector notation V = 
(do, A x , A y ) etc., with lattice derivatives in the space di- 
rection. Note, there is no charge Chern Simons term since 
the net Chern number vanishes. The first term arises be- 
cause of the opposite Chern numbers for the two spin 
bands. It associates a spin of K/2 with the 2-7T flux of 
a. Finally, the dynamical gauge field (a) is non-compact, 
i.e. there are no monopole configurations created by the 
dynamics, since a monopole event (which inserts 27r flux) 
is associated with a spin flip, which is forbidden by S z 
conservation. We now integrate out the dynamical gauge 
field a to obtain: 

Ce[A1] = Z&K [Al? + f [(V X ~2 AS)2] (2?) 

the first term is a London term that indicates the phase 
is a 'spin superfluid'. This implies XY magnetic order, 
although it does not specify the precise pattern of order- 
ing (i.e. ferromagnetic, vs various Neel ordered states). 
However, given that the mean field theory is favored by 
the microscopic Hamiltonian Eq. (|2"D|) , it is natural to 
expect this to be the usual Neel ordered state with op- 
posite moments on the two sublattices. Note, since we 
have already doubled the unit cell, the Neel order is at 
q = (0,0), but transforms nontrivially under the point 
group. 

The nature of the order may be established more rig- 
orously as follows. Note, due to the absence of monopole 
events in the dynamics, one is left with a Maxwell ac- 
tion for the dynamical gauge field a. The gapless pho- 
ton excitations that this implies, corresponds to oscilla- 
tions of the dynamical magnetic flux. However, due to 
the binding of flux to spin density, this implies a fluc- 
tuating S z density. In fact, this gapless photon is sim- 
ply the Goldstone mode associated with XY spin sym- 
metry breaking (see the discussions in section IV). In 
order to establish the spin order one needs to evaluate 
(S+Sq). This corresponds to evaluating the correlators 
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of the monopole insertion operator (V^Vb). It is well 
known^i, that in magnetically Coulomb phase in D=2, 
the monopole operators have long range order. The pre- 
cise ordering pattern corresponds to the transformation 
properties of the monopole insertion operator under the 
lattice symmetries. This is explicitly evaluated in Ap- 
pendix [F] where it is confirmed that we indeed obtain the 
Neel state. Note, this phase is entirely identical to the 
regular XY Neel ordered state, although it is obtained 
from such an 'exotic' starting point. The magnetic order 
induced by gauge fluctuations ensures, for example, that 
the counter-propagating edge states implied by the mean 
field theory are gapped. 



B. Charge Fluctuations 

Consider lowering the charge gap to approach the Mott 
transition so that the electron occupation at each site 
can now fluctuate into the n = and n = 2 state as 
well. This is achieved by lowering the value of U in equa- 
tion 1191 These fluctuations are typically incorporated 
via the slave rotor formulatio n 30 ' 31 , by introducing a ro- 
tor variable z r = e l ^ r and its conjugate number variable 
N r at each site. The electron operator is then written as 
c\ a = z r fl a) and the constraint now reads n( + N r = 1. 
This however does not capture all Mott transitions, for 
example it was pointed out that Hubbard models with 
psuedospin symmetry requires introducing a SU(2) slave 
rotor— which is a simple extension of the SU(2) slave 
boson theory 3 - 3 * to include both electron and hole fluctu- 
ations. A pair of complex fields introduced, 
which are SU(2) 'rotor' variables i.e. \z lr \ 2 + \z 2r \ 2 = 1. 
The electron operator is written as: 



c r] — frf z rl Vrfli z t2 
c rl = frl z rl + ?7r .f^T 2 r2 • 



if we define: 



Z r 

9r 



Zil -z 
Zi2 4 



Fr 



Vrf^ 



(28) 

(29) 

(30) 
(31) 



the electron operator may be written compactly as in the 
last line; where Z r is an SU(2) matrix. The mean field 
Hamiltonian (|25[) takes on a particularly simple form in 
this notation: 

H MF = F r ' + l rr' + tfrT>l*z)Ff (32) 



with Xrr' = X^r'+ixlr' an d with diagonal spin dependent 
hoppings t s rr , as shown. 

This definition has an SU(2) redundancy, clearly any 
set of SU(2) matrices U r can generate the transforma- 
tion Z r — » Z r Ul, F r — > U r F r leaving the form above 



invariant. Naturally, there is a close connection between 
this redundancy and the constraint that needs to be im- 
plemented to obtain the physical Hilbert space. If the 
operators T r generates the above transformations, then 
the constraint on physical states is that T r = 0. 

Luckily for our purposes we will not need the full 
SU(2) formulation. Instead, we only retain the fact 
that there are a doublet of charged bosons, but as- 
sume that they are rotor variables with equal amplitude 
Z\ f = -^ e ^ lr ') z 2r = ~^2 e ^ 2r - This form is necessary, 
and sufficient, to incorporate the gauge fluctuations and 
symmetries of the mean field theory Eq. (f25|) . The SU(2) 
matrices U that preserve this structure are generated by 
U e — e 16 ^* , and U' — ifj, x where \i are 2x2 Pauli matrices 
in the usual representation. While the former is simply 
the 17(1) gauge transformation that leaves the mean field 
Hamiltonian invariant, the latter appears as the gauge 
transformation required to implement reflection symme- 
try (Appendix [C]) . Hence, in the absence of a two com- 
ponent slave rotor theory, we would be unable to im- 
plement this physical symmetry starting with the mean 
field theory Eq. ([2!>]) . If we now define N\ r and N 2r as 
the operators conjugate to the two phase variables, the 
physical states in the basis |JVi, N2; nt, n\ > satisfy the 

constraints (Ni + AT 2 ) + (nt +n\) = 1 and invariance un- 
der the combination of |JVi, N 2 >-> (-l) Nl \-N 2 ,-N 1 > 



constrains 



and \n f >^ r) r (f r] f rl + /J ; /L)|n/ >. This 
the physical states on a site to be: 

I T> = |0,0;T> (33) 

I l> = |0,0;j> (34) 

|0> = (|0,1; >+7j r |-l,0;U>) (35) 

I U> = (-»?r|l,0; >+|0,-l;U>) (36) 

A physical interpretation of the charge boson doublet is 
provided at the end of appendix IG 31 

The superconducting phase 

An advantage of the two component slave rotor formu- 
lation above is that the superconducting state is very nat- 
urally described. As the Mott transition is approached, 
the gap to the charge excitations closes, at which point 
Z\ , zi condense. If both bosons condense individually, 
< z\ >— $1, < 22 >= < I > 2 a superconductor is obtained. 
This can be verified by noting that the gauge invariant 
quantity < z\z 2 >= &1&2 which is the superconducting 
order parameter, has long range order. 

Once the charged bosons are condensed, the electron 
^ r has a finite overlap with the 'spinon' F rcr , and in- 
herits its dispersion. Since ^> r ~ (Z)F r , and (Z) = 



\®2\ 2 Uq,, where {7$ is a unitary matrix, an ef- 



fective Hamiltonian that describes the electron excitation 
can be readily obtained from the mean field Hamiltonian 
El: 



(37) 
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when both bosons are condensed with equal magnitude 
$! = e l6, $ 2 , the electronic Hamiltonian is just: 

HdSc = ^2 C i<r{Xrr> + Kr'Yr'o' + A^e^/C^C^, + h.C. 
rr' 

(38) 

where the pairing function A rr / = ie l6 ri r >xi r i is sym- 
metric upon r <-> r' and changes sign under ir/2 ro- 
tation. Hence it is a spin singlet d-wave pairing func- 
tion. If t s = this is the regular d x 2_ y 2 superconductor 
with translation symmetry and nodes at (±7r/2, ±7t/2). 
The addition of the spin dependent hopping opens a 
gap at these nodes. It doubles the unit cell and breaks 
time reversal symmetry (although the combination of 
a unit translation and time reversal remains a symme- 
try). Note however, since the above SC Hamiltonian 
is a unitary transformation of the hopping Hamiltonian 
Hmf, it shares the same spectrum. In particular, since 
the hopping Hamiltonian had a nontrivial band struc- 
ture and was argued to have counter-propagating spin 
filtered edge modes. The same is true of the supercon- 
ducting Hamiltonian above. However, once the charge 
bosons condense, additional terms are allowed in the 
fcrmion Hamiltonian, which are forbidden in the insulat- 
ing state. For example, one can add terms corresponding 
to 'spinon pairing', which are forbidden by the gauge 
U(l) symmetry on the insulating state. Such terms can 
arise simply from the electron's next neighbor hopping, 
which in the two component notation corresponds to 
ipl^ipr' = F ] r ,Z^ z ZF r = \<P\ 2 Flfi x F r . Note however, 
the strength of such a term is proportional to the conden- 
sate density, and becomes very small near the transition. 
Including these terms removes the nontrivial topological 
properties (such as gapless edge states) of the supercon- 
ductor. Therefore, the superconducting state obtained 
on boson condensation is a generic d-wave superconduc- 
tor with doubled unit cell. If, however, we demand par- 
ticle hole symmetry, then terms like the second neighbor 
hopping term above are forbidden, and the superconduc- 
tor indeed has nontrivial band topology (see appendix 
IG 2[) . In our previous analysis it was convenient to con- 
sider a particle hole symmetric superconductor, although 
it is difficult to demand that as the microscopic symmetry 
of a physical system. Luckily, we show in Sec lVIII that de- 
manding complete particle hole symmetry is unnecessary 
- instead we need only tune one chemical potential type 
term to zero, to obtain a direct relativistic transition. 

We note in passing that condensing just one of the two 
(z%, Zz) fields leads to a topological band insulator (TBI). 
(At the end of Sec IV CI we show that these states are real- 
ized in the SC vortex core.) In choosing one of these fields 
to condense, one spontaneously breaks the reflection sym- 
metry. If z\ condenses, the appropriate mean-field elec- 
tronic Hamiltonian can be obtained from Eq. (|25p by 
replacing f a — > c a . Similarly if z 2 condenses we simply 
make a different replacement: f a — > r] r e arT i c^, . The above 
different condensation scenario of z\ and Z2 can be sum- 
marized by the vector ft = z^az where z = ({ l Zi), (^2))- 



1 DL j 














SC / 



TBI 2 

FIG. 6: The ft sphere and the corresponding electronic phases. 

Using ft we represent the topological insulator (TBI) as- 
sociated with (z\) ^ and (z 2 ) = (TBIi) by ft pointing 
to the north pole. Similarly the TBI associated with 
(zi) = and (z 2 ) ^ (TBI 2 ) by ft pointing to the 
south pole. The particle hole transformation transform 
TBIi into TBI2 and vice versa. The d-wave SC discussed 
throughout the early part of the paper is represented by 
ft lying in the xy plane. This is illustrated in Fig.f6|). 

C. The field theory of transition 

Consider approaching the transition from the insu- 
lating state. The charged bosons are gapped but ac- 
quire a dispersion from the electron hopping Hamiltonian 
[T9l This includes the regular nearest neighbor hopping 
■^hop = — ^S(rr') c ra c r'cr'7 which we treat in detail below 
(other hopping terms lead to similar results). Expand- 
ing the electron operator as in Eq. (|28p. and assigning 
the mean field expectation values to the spinon bilinears 
J(12<j frafr'v) = Xrr' , one also obtains a staggered flux 
dispersion for the z-bosons: -f/hop = — j S Xr'r(z\ r Zlr' + 
z\ r Z2r')- This dispersion has a unique minimum (for flux 
$ 7^ 7r) which is relevant to the low energy theory. 

The low energy field theory of the transition is argued 
below to take the following form: 

C E = C z + C f (39) 

Cf = i—A b ■ V x a 
1 2tt 

V = r(|z 1 | 2 + |z 2 | 2 )+C/(| Zl |2 + | Z2 |2)2_ A | Zi |2| Z2 |2 

AC = ix[z*{d T — iaa)z-L — z^ipr — iao)z-i] 

The Euclidean Lagrangian consists of two pieces, one 
arising from integrating out the fermionic spinons C f and 
one for the z-bosons C z , coupled to an low-energy U(l) 
fluctuating gauge field a. Note, the representation of the 
electron operator in Eq. [ITI] implies that the z-bosons 
are also minimally coupled to the U(l) gauge field of the 
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spinons. The remaining terms appearing in C z are the 
ones allowed by the symmetries of the problem, which 
act on these bosons in a nontrivial manner as described 
in Appendix [C] Since the fermionic spinons are gapped 
throughout the transition, they can be safely integrated 
out. This produces the spin Hall effective action Cf 
since the mean-field state of the spinon is a spin hall 
insulator. (Note that to the spinons a plays the role 
of the charge gauge field.) Symmetries prohibit other 
terms in this low energy action. For example, (see ap- 
pendix[C]) under reflections about the x or y axis passing 
through a plaquette center, the fields transform accord- 
ing to Z\ — > z|> z 2 — * — z*, in addition to the reflec- 
tion of coordinates. This symmetry rules out terms like 
\ z i 1 2 — \ z 2 1 2 and z{d T zi + z\d T Z2. It does however al- 
low the linear time derivative term AC. While point 
group symmetries do not forbid a background flux term 
like h(d x a y — d y a x ), taking into account symmetry under 
time reversal (followed by translation) rules out such a 
term. (Only a staggered background flux is allowed, but 
this has already been incorporated while deriving the z- 
boson dispersion.) 

In the absence of AC, the effective action for the tran- 
sition is in the easy-plane NCCPi universality class^. As 
we discussed at the beginning of the paper, the nature 
of direct transition (whether it is generically fluctuation 
driven first order) in the relativistic easy-plane NCCPi 
model is controversial. In contrast, in the presence of AC 
the non-relativistic NCCPi transition in d + z = 4 has 
been studied using perturbative RG 34 . The conclusion is 
that fluctuation generically drives the transition to first 
order. 

In order to access the relativistic transition, one needs 
to tune both r and \x to zero in equation (|39[l . The lat- 
ter corresponds to a chemical potential since z\ and z^ 
carry opposite electrical charge. It is important that one 
only need to tune two parameters to reach the relativis- 
tic critical point. Because it implies this critical point is 
realizable in a two dimensional phase diagram. Physi- 
cally jj, measures the deviation from half-filling. In the 
AF insulator phase the charge density remains fixed for 
a range of fx, due to the charge incompressibility. In the 
SC phase the charge density does change with chemical 
potential and one would need to fine tune fi to reach the 
relativistic NCCPi quantum critical point. 

The above situation is analogous to the simpler case of 
accessing the relativistic XY transition in the superfluid- 
Mott insulator transition of the Bose Hubbard model 35 . 
There, the relativistic transition occurs at the tip of the 
superfluid lobe, while moving away from this point leads 
to a non-relativistic Bose Einstein condensation transi- 
tion with dynamical critical exponent z — 2. In Fig. ([7]) 
we present an analogous (schematic) phase diagram in 
the plane spanned by the chemical potential /i and inter- 
action strength U . 

As being discussed in Sec.V A, the electromagnetic 
field associated with the gauge field a corresponds to 
the 3-current of S z arising from the gaussian fluctua- 



tion in the magnetic XY order. Since z% and z 2 couple 
to a as charges, it is natural to associated them with 
the magnetic vortices. From the earlier discussions we 
see that they carry electrical charge. A related observa- 
tion is that in the superconducting state near the critical 
point, the superconducting vortices carry spin. This can 
be seen as follows. Consider an equal amplitude conden- 
sate (zi) — y/pe 1 ^ 1 , (z 2 ) = y/pe 1 ^ 2 . The energy cost of a 
static vortex is then obtained from the lagrangian above 
as£ = /o(V0i - a - eA/h) 2 + p(V(f> 2 - a + eA/h) 2 , where 
A is the external vector potential that induces vorticity. 
A unit vortex with § A ■ dl = h/2e can be created either 
with (1) § V0i • dl = 2ir, § V</>2 • dl — and § a-dl = w or 
with (2) § V<j&i -dl = 0, § V0 2 • dl = 2tt and § a ■ dl = -n. 
Due to the presence of Cf which attaches spin to emer- 
gent flux, these vortices carry spin ±ft/2. Thus, near the 
transition, the low energy vortices in the antiferromag- 
net carry charge ±e while vortices in the superconductor 
carry spin ±h/2. Thus, a direct transition between these 
phases can be explained on the basis of the condensation 
of these defects. 

Finally before closing this section we note that the 
mechanism that generates magnetic order studied here 
does not require "confinement" which is invoked in other 
slave boson theories. The advantage of our mechanism is 
that it allows us to identify both the nature of the mag- 
netic order, which has not been possible in confinement 
type approach, and the superconducting pairing symme- 
try unambiguously. 



VI. ROAD MAP TO THE APPENDICES 

Before concluding, we outline the content of a number 
of appendices. They are presented either to elaborate 
points made in the text, or because they are interesting 
digression which contribute to the general understanding 
of subjects covered by this paper. Appendix |A1 contains 
the technical details of the duality transformation. Such 
transformation allows one to go back and forth between 
Eq. © and Eq. (filjj) . Appendix [B] contains the Bogoli- 
ubov deGennes treatment of the core state of supercon- 
ducting vortices. It compliments our numerical results 
presented earlier. Appendix [Dl present a mean-field the- 
ory for the XY ordered state and study the core states 
of the AF vortices. Appendix lEl provides details of the 
vortex tunneling amplitude calculation in both the SC 
and the AF phases, where we explicitly compute the XY 
order pattern (AF) from the SC vortex tunneling, and 
the SC pairing symmetry (d-wave) from the AF vortex 
tunneling. Appendix [C] is devoted to the symmetry anal- 
ysis of the slave boson theory in the AF phase and the 
transformation properties of the antiferromagnetic vor- 
tices. Appendix IPl presents the quantum number of the 
monopole operator of Eq. (Ti"6|) within the slave rotor the- 
ory. Appendix [Gl gathers some miscellaneous discussions 
mentioned in the main text. Finally Appendix [H] dis- 
cusses the direct transition between an f-wave SC and 
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AF magnetic phase on the honeycomb lattice, which is 
a direct application of the techniques that we learned in 
the current square lattice case. 



VII. CONCLUSION 

In this paper we present a theory exhibiting a d-wave 
superconductor and an easy-plane antiferromagnetic in- 
sulator as its two phases. Unlike ordinary d-wave su- 
perconductor, our superconductor has a full gap due to 
nodes mixing caused by a translation symmetry break- 
ing. In addition, the topological defects, namely vor- 
tices, of each of these phases carry the quantum number 
of the order parameter for the other phase. We derived 
the effective field theory describing the superconductor 
to AF insulator phase transition. It is the easy-plane 
non-compact CPi model. The possible phase transition 
scenarios of this model, include: (1) a direct continuous 
SC-AF transition, (2) a direct first order SC-AF phase 
transition, and (3) two consecutive continuous transitions 
linking SC to a SC+AF coexistence phase and finally 
to AF. In particular, if scenario (1) is realized, it would 
constitutes an realization of Landau forbidden transition. 
While the first two scenarios correspond to the conden- 
sation of the fundamental vortices in theory, the third 
requires the condensation of a composite of the funda- 
mental vortices. 

Physically this phase transition requires tuning the in- 
teraction strength (or the bandwidth) at fixed, half, fill- 
ing factor. We specified the symmetry requirements for 
its existence. When the filling factor deviates from 1/2 
the transition generically becomes first order. In Fig[7] 
we present a schematic phase diagram in the interac- 
tion strength - chemical potential plane. For a range 
of chemical potential where the superconductor remains 
fully gapped with the minimum gap (left panel of FigfS]) 
occurring at (±7r/2, ±7r/2) in the square lattice Brillouin 
zone. For sufficiently large chemical potential deviation 
(from that of half filling) the superconducting gap closes. 
When that happens the quasiparticle spectrum becomes 
that of a conventional d x 2_ y 2 superconductor (right panel 
of FigE]), and the positions of the node change with chem- 
ical potential. 

In the above discussion we have assumed that upon 
doping the AF-SC transition is direct. It is also possible 
that such transition proceeds in an indirect manner. For 
example, imagine a situation the doped charge in the AF 
insulator are accommodated as charged vortices zi and 
Z2 ■ To minimize the cost in kinetic energy, it is favorable 
to keep the total vorticity zero, so that there are equal 
number of vortex and antivortex. Thus, if the doping is 
n-type, there will be an equal number of z\ and z\ vor- 
tices (recall that Z\ and Z2 vortex carry opposite charge), 
while for p-type doping there will be the same number 
of 22 and z\ vortices. Due to the logarithmic attraction, 
it is energetically favorable for z\ , z\ or z\ , Z2 to form 
charge ±2e bound pairs. Because these pairs carry no 
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FIG. 7: The schematic global phase diagram. The horizon- 
tal axis is the interaction strength, U, and the vertical axis 
is the chemical potential /i. The dashed lines mark the non- 
relativistic phase transition at mean-field level (which is likely 
to be first-order beyond mean-field theory). This transition 
line terminates at the PH symmetric critical point P* dis- 
cussed in the present work. Along the dotted line the average 
charge density is (n) = 1, i.e., half-filling. The SC phase 
surround the the AF phase is a fully gapped d-wave super- 
conductor. Upon changing of chemical potential the quasipar- 
ticle gap closes and the SC undergoes a second-order Lifshitz 
transition into the nodal gapless d-wave SC at the dash-dotted 
lines. The quasiparticle spectrum of the gapless d-wave super- 
conductor is identical to that seen in high T c superconductors. 




FIG. 8: In the presence of a sufficiently large chemical po- 
tential fj, the nodes of a conventional d-wave superconduc- 
tor is restored. The parameters used in this figure are 
4> = 7r/5, t = 0.2 and left panel fi — 0, right panel /i = 0.5. 



net vorticity, their condensation does not destroy the AF 
order. Hence a SC and AF coexistence phase emerge. 
Such a transition is of the usual non-relativistic XY uni- 
versality. In the coexistence phase, due to the non-zero 
{z\z\) or (z\z2), the z\ and Z2 AF vortices mixes, and 
as the result only one AF vortex z — az\ + (5z2 (where 
a, [3 are complex mixing coefficient) is left as low energy 
excitations. When z condenses the AF order is finally 
destroyed and the system becomes a pure superconduc- 
tor. The coexistence — * AF transition is also of the usual 
XY universality class. 

It is also interesting to look at the indirect transition 
from the superconducting side. As discussed at the be- 
ginning, such transition is triggered by the condensation 
of the composite vortex ^/jst^. Such composite vortex 
carries spin quantum number 1 but no net vorticity (of 
the superconducting order parameter). As the SC to 
SC+AF coexistence transition is approached, this excita- 
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tion will become the lowest-energy magnetic excitation. 
Like the roton in a superfluid, such vortex-antivortex pair 
will carry a finite momentum. Such momentum will shift 
the momentum of the minimum energy ^jvf'j. away from 
(tt, 7r), similar to the incommensurate magnetic excita- 
tion of the high T c compounds. 

We note that the similar ideas that the defects in the 
AF phase can carry charge and the defects in the SC 
phase can carry spin, are introduced in the phase string 
theor y 10 ! 11 ' 12 in the context of High-Tc superconductor. 
Although there are some significant differences between 
the current work and the phase string theory: for ex- 
ample, we require the doubling of unit cell, easy-plane 
spin anisotropy and full energy gap on both AF and SC 
phases, it is still possible that there are underlying rela- 
tions which are left for future investigations. 

The above discussions are intended to point out some 
similarity of the phenomenology in our model and that 
of the cuprate superconductors. However, due to the 
symmetry difference in the superconducting state, the 
connection is by no means clear. Nonetheless it is inter- 
esting to explore the possible relation with the cuprates 
in future studies. 

This work is supported by NSF-DMR 0645691, LBNL 
DOE-504108 and DOE DE-AC02-05CH11231. 



APPENDIX A: DUALITY AND SYMMETRIES 
OF THE SUPERCONDUCTOR VORTEX 
THEORY 

We first fix the relation between the relativistic vor- 
tex fields VPj, and the canonical vortex creation and 
destruction operators. In Eq. 

*T = (^AV T + e^BVi 



2 
2. 



(Al) 



where i/av t bosonic vortex field for the AVf 



vortices, and 

etc (For more discussions see the end of this paragraph) . 
The phase (3 is the phase of the annihilation matrix ele- 
ments of (AVt, BV |) and (BVj, AV T), respectively (the 



reason that these two phase factor are the same is due to 
the To t^l symmetry). We note that in terms of Wf.x 
Eq. (|9|) has a relativistic form. It can be derived from 
a theory where ^AVf , V'AVj > V'bv^ and 4'bv 1 are treated 
as non-relativistic fields; in the presence of pair cre- 
ation and pair annihilation terms — | J|e l ' 3 '!/'AV T '0BVx + c - c 
the combination (i/'avt — e_ *'vBVl)A/2 and (V'bvi — 

e_I/3 V'avt ) become the more massive fields and hence 
can be dropped from the low energy theory. The left 
low energy combinations is given in Eq. (jAip . The im- 
portant thing is that in the derivation sketched above it 
is assumed that Vav t > VaVj > '<Pbv 1 and ipBV l are a U Re- 
generate with respect to one another (in the sense that 
permuting any two fields will leave the (non-relativistic) 
action invariant. This important degeneracy is guaran- 
teed by the symmetry of operations of Eq. ([5]) . 

Note that according to Eq. ([5]) PH transforms a vortex 
into an anti-vortex, while preserving spin, while mirror 
reflection V also transforms a vortex into an anti-vortex, 
while preserving the spin. Finally, T o TTZ transforms a 
vortex into an anti-vortex and flips spin. Thus we have: 



PH: AVf 
=^* T - 
V: AVf 

T o TTZ: AVf 
- 



AV 

t 



I 

AVf 



t 



BV 

t 



T 

BV X , 

t 

7>(*)T 

AVf, 
t 



(A2) 
(A3) 



\&T (anti-unitary) 



(A4) 



this allows us to restrict our field theory to the form in 
equation Eq. ©. 

Duality: 

We now carry through the duality starting with 
Eq. ©. We first set ^ a = \fp4>„ at long wave- 
length/time, where p is the average condensate density of 
the up and down vortices, and is a U(l) phase factor. 
Subsequently we can introduce auxiliary fields J° so that 



J 



J I'tJ;]^"^ 2 *^^^ 2 "^^^^^^^ (A5) 



r 



Next,we separate 4> a d^(j)a into 4>ad^4> a = id^Oa+^d^, 
where <j) v a denotes the topological non-trivial part (i.e., 
the vortex containing part) of <j) a . Integrating out Re- 
generates the constraint d^J® = 0. This constraint is 



solved via the introduction of two new "gauge fields" = 
^(Vxa"),,. Now we collect all the and a M dependent 
terms in the action to obtain 
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1 

2~7T 



-(Vxa) r -(Vx4 



-(V x of +i-A- (V x a) }. (A6) 

2 7T I 



In Eq. jMJ) K° = -^Id^l is the 3-current of the 
vortex of superconducting vortices (i.e., the vortex of ip a 
bosons). Next, we integrate out a M which leads to 

/J(V x a°") M = 2e(V x A)^ + less relevant terms 

a 

at long wavelength. (A7) 



Eq. (|A7[) hxes the flux in + a* and leaves the flux of 
a 1 — free to fluctuate. Let us define a M = (a^ — ajJ/2, 
then the remaining part of Eq. (IA6|) reduces to 



Sd 



ual 



( d 2 xdt^2 j^[e 2 (V x A)£ + (V x a) J] -ieA^JfT + +ia„(K^ 

mi. ** ' 



(A8) 




FIG. 9: The reduced Brillouin Zone (shaded area) of the 
model Eq.([T}. 



If we identify 

the 3 — current of Z\ boson 
the 3 — current of boson, (A9) 

Eq. (|A8|) can be recognized as the Feynman path integral 

I 



representation of the second quantized action in Eq. (IT 



APPENDIX B: BOGOLIUBOV DEGENNES 
ANALYSIS OF THE SUPERCONDUCTING 
VORTICES 

In section U we have determined the core structure of 
the superconducting vortices numerically. To make sure 
that our numerical results are universal, in this section we 
present an analytical study in the long wavelength limit. 
Interestingly, as a by-product, we obtain a new situation 
where zero modes exist in a non-Dirac like Hamiltonian. 

We first write down the d-wave pairing Hamiltonian 
Eq.(TTJ in the momentum space: 



(Bl) 



where 



/ 2t(cos fe 2 - cos k ± ) xi(l + < 

XI (1 + e*< fc l+ fc 2) + c ifc l + c ik 2) 

X2<-1- 



i(fc 1 +fc 2 ) + e _ife l + c 
-2t(cOS k 2 - cos fcj) 
-»(fel+Jt 2 ) + e- ifcl + , 




k 2) X2< _1 ~ e- i < fc l+'=2) + , 

ik 2) 2t(cos* 2 - cosfei) -xi(l + e - i < fc l+ fc 2) + t 

-Xl(l + e i < fc l+ fc 2) 4. e «fel + e «*2) -2t(cos fc 2 - c < 



i hi) 



*2)\ 
k 2 ) 



(B2) 



r 



In Eq. (JBTj) and Eq. (JB2J) 



and we denote Re{\) — Xi an d I m (x) — X2- We allow xi 
to be a complex number since the pairing amplitude gain 
(B3) a phase under the charge-£/(l) rotation. k\ and &2 are 



15 



the momentum components along the reduced reciprocal 
lattice vectors in Fig.@. 

Let us first study the parameter regime: ~ 1, 

|X2 1 ~ 1 and \t\ « 1. If t = we have a gapless d-wave 
superconductor with two inequivalent Dirac nodes at K 
and K* (see FiglSJ). The t term mixes K and K* and 
opens up a small gap. In this limit we can linearize the 
Hamiltonian matrix Eq. (|B2[) around K as: 




At 






-At 
2xiik 2 




where (fci, k 2 ) = k — K. 

In the presence of a vortex \i 
Eq. (|B4p . in real space, becomes 



X2(x) = \x2\e~ 



(B4) 



id(x) 



H 



K 



-Atfi 3 + 2ixiH2T 3 d 2 - 2i|x 2 |e l2T3 M 2 Tie 8 2 T3 di. 



(B5) 



Here the Pauli matrices \x mix the sublattice labels and 
r mixes particle and hole. One can perform a series of 
rotations to transform Hk into a more familiar form (a 
form that has been studied before^): 

H K = At(cos9r 1 + sin0r 2 ) - 2i(|x2|mit 3 <9i + XiA i 2T 3 9 2 ) 
=4i(cos0Ti + sin 0t 2 ) - 2«(^it 3 c\ + fi 2 T 3 d 2 ) (B6) 

where we define d\ = 2|x 2 |9i and d 2 = 2xid\. 

The transformations that lead from Eq. (|B5[) to 
Eq. (|B6|) are the following. Firstly we perform a rotation 
along t\ such that r 3 — ► — r 2 . Afterwards we perform 
another rotation along \i\ such that /i 3 — > — /i 2 , /i 2 — > /i 3 . 
After these transformations 



-> 4t/i 2 - 2ixiT 2 fi 3 di - 2i\x 2 \e 



-i9/ 2 T 2 



T\[i 3 e 



(B7) 



-i6/2T2 e i8/2 l i 3 _ Al- 



Second, we perform the rotation e 
though this rotation is spacial dependent, it commutes 
with d in the long wavelength limit. Another worth not- 
ing fact is that e - l6 / 2T 2 e ^/2fi 3 nag n0 branch cu t since 
for both = and 9 = 2tt, e -^/2r 2e te/ 2f , 3 = L After the 
above transformation 



Hk — ► At(cos9/j, 2 - smO/ii) - 2i(xir 2 // 3 9i + 



X2klM33 2 ) 

(B8) 



After interchanging \i «-> r and redefining — > 9 + n/2 
we obtain Eq. (|B6p . 

In Ref. 3S it is shown that Eq. (|B6j has one zero mode. 
Similarly one can show there is another zero mode as- 
sociated with K* . And same analysis can be applied to 
anti-vortex. We thus established the two mid-gap modes 
in the SC vortex (anti- vortex) . Their existence ensures 
the existence of the four low energy spinful vortices dis- 
cussed in Sec|U 



Next, we analyze a different limit of parameters: |xi | ~ 
1, \t\ ~ 1 and |x2 1 "C 1. Since it is possible to adiabati- 
cally tune the parameter from the previous regime to this 
limit, while preserving 90° rotation and PH symmetries, 
we expect the zero modes to survive. In specific, in the 
previous limit the two zero modes are associated with 
two different nodes (for small t and slowly varying #(x) 
momentum is approximately conserved), and transform 
into each other under 90° rotation. By forming sym- 
metric and anti-symmetric combination of the two zero 
modes we obtain two new ones that transform under 90° 
with s and d symmetries respectively. On the other hand 
PH does not change angular momentum, and transforms 
each mode into its particle-hole conjugate hence reverse 
the sign of energy. Since the entire adiabatic process pre- 
serves the 90° rotation and PH symmetries, the in-gap 
states can not shift away from zero energy. We have per- 
formed numerical study of the vortices in this parameter 
limit and find the mid-gap modes. 

The reason we are interested in the second limit is be- 
cause it gives rise to an interesting situation where zero 
mode arise from a non-Dirac like cone. In the limit of 
|X2 1 - * two quadratically dispersed bands touch at a 
single X point (see Fig[9]). Nonzero |xa| opens up a gap. 
To quadratic order in momentum departure from X the 
Hamiltonian looks like: 



H 



x 



(t(k\-k\) 
-X\k\k 2 


-4X2 



-Xik\k 2 
-t(k 2 — fcj) 
-4X2 






-4x 2 
t(k 2 — fci) 
X\k\k 2 



-±x* 2 



X\k\k 2 
-t(k 2 — k\)/ 
(B9) 



where (fci, fe) = k~X . The existence of zero modes when 
the is a vortex in X2 is very interesting. For up to present 
all known zero modes are associated with "mass vortex" 
in Dirac-like equations. The fact that Eq. (|B9[) is 11011- 
relativistic, yet in the presence of "gap vortices" there are 
zero modes suggest a new type of "index theorem" : For 
one vortex in xi of Eq. (|B9[) there are two zero modes. 



APPENDIX C: SYMMETRY ANALYSIS FOR 
THE SLAVE ROTOR THEORY 

In this appendix we study the symmetries of the slave 
rotor theory. As discussed in appendix [Vj the relation 
between the electron operator and the slave rotor and 
spinon operator is given by^i: 





(CI) 



or equivalently Eq. ([28 



(C2) 
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According to the above equation one can insert a site- 
dependent SU(2) matrix and its inverse between the /- 
spinon and z-rotor matrices while leave the electron op- 
erators invariant. In terms of /-spinon and z-boson this 
local SU(2) transformation is: 



«•; = ( ^ 

VHJil 



Zi = 



Zil 



i«i.§ . z 



(C3) 



ifestation of various physical symmetry of the AF phase 
in the slave-rotor gauge theory. 

There are two physical global U(l) symmetries: the 
global spin (S z ) rotation, and the charge-[/(l) symmetry. 
They are manifested in the slave-rotor theory as: 



S z rotation by 8: fa 



fs 



-iO/2 f 
Ji. 



i(9/2 



Cil, 



(C5) 



The above SU(2) gauge symmetry is broken by the form 
of the mean-field spin-Hall insulator Hamiltonian of the 
spinon (Eq. 10). The remanent gauge symmetry is U(l), 



and 



fia 



7i, 



(C4) 



Due to this remanent U(l) symmetry, the fluctuation 
around the mean-field theory appear in the form of a U(l) 
gauge theory. We have shown in appendix [V] that if the 
U(l) gauge field is in the Coulomb phase, the spin-Hall 
insulating properties of the spinon indicates that this is 
actually the XY ordered phase. This is because the spin- 
Hall response implies that the photon of the gauge field 
is actually the Goldstone mode of the XY ordered, as has 
been shown in text. In the following we study the man- 



■ ) ■ fia 


~* fia 








iO 




-10 


Zil - 
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Zi2 ~ 


■+ e z, 2 




i6 




id 


=>■ Cii - 


-» e Cjf 


Cil - 


■* e di . 



(C6) 



These equation imply Z\ and Z2 are spin zero charge- 
carrying particles. 

Next we study discrete symmetries. Before doing so it 
is convenient to write the mean-field spinon Hamiltonian 
in the SU(2) form. In this form any transformation which 
changes the mean-field Hamiltonian by a SU(2) gauge 
transformation is regarded as a symmetry operation due 
to the gauge redundancy. Eq.Q can be rewritten as: 



J 



-x*j o 

~Xij 



«*J'» 



o -u 



r 



ipj + h.c. = iplUiji/jj + h.c 



(C7) 



Under a SU{2) gauge transformation ,i.e.,ipi — > Witpi 



U, 



W]Ui 3 W r 



The SU (2) gauge invariant quantities are 

Tr [U Xi , i+t Ut i+ ,, i+ ,U Xi+ii% ] = -2icos20. 



(C8) 



(C9) 



for every triangular loop i^i + x—>i + y—>-i. Two 
mean-field Hamiltonians with the same loop trace are 
related to each other by an SU(2) gauge transformation. 

Now we are ready analyze the discrete symmetries. 
First, it is obvious that the 90° rotation around the cen- 
ter of each plaquette TZg® is a symmetry since it leaves 
Eq. (0 invariant. Under this operation the z bosons 
transform as: 



n 



90 : Zn — > Z7^ 90 (i) ! i 



Zil 



Zn 90 (i),2- 



(CIO) 



Next we consider the reflection V about the verti- 
cal/horizontal lines passing through the center of each 
plaquette. Naively one would claim that this is not a 
symmetry of Eq. @ because 



But this is actually a symmetry transformation because 
it does not change loop trace; one just needs to per- 
form a global SU (2) gauge transformation Wi = ia% (see 
Eq. ([U3)) ) to restore Eq. ©. Thus reflection, V, is a sym- 
metry. Under V the z bosons transform as (see Eq.JC3|): 



WoV: zn 



V(i),V 



(C12) 



This is consistent with the result Eq. (|G2p . Next, come 
the particle-hole transformation PH. According to 
Refill, this transformation is implemented on the z- 
boson via 



PH. jioi > fia 
^>Cil 



m4i 



Zil 

Cil 



Zi2 



Zi2 



-f}i c il 



(C13) 



Again, this transformation law is consistent with result 
in Eq. (|Gl[) . This transformation also preserves a gapped 
z boson action and Eq. @ hence is a symmetry. Now we 
consider the translation T: 



V: X 



X 



t -> t . 



(Cll) 



T-xfy fi * fi-\-xJyi %i ^ %i-\-x/yi ^icr * Ci-\-xfycr- (C14) 
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Under this operation the parameters in Eq. ([2]) transform 
as: 

T: X ^X\ (C15) 

This causes the loop trace to change sign (i.e., 
-2icos20 -> 2t cos 2(f)) hence T is broken by Eq. ©. 
Next, we proceed to time reversal TTZ. Under time re- 
versal 

T7?.(anti-unitary): f ia -> {io- 2 ) a pf il3 , z ia -> z ia 

=^ Q T ^C 4i , Cj; -Cj| (C16) 

Under the above transformation the parameters of 
Eq. {2]| change according to 

Tft: X ^X*, (C17) 

which also flips sign of the loop trace. Hence T72. is also 
broken by Eq. $Z§. 

Next we study the compound operation ToTIZ . Com- 
bining Eq. (|C15|C17|) it is simple to show that Eq. ^ is 
invariant under this transformation. Thus we conclude 
T o TTZ is a good symmetry. Using Eq. (|C14IC16|) we 
obtain the z boson transformation law as 

To TTZ : zn -> z l2 -> Zi+i,2- (C18) 

Once again, this is consistent with Eq. (|G3[) . 

Because the slave-rotor theory does not capture the 
S z rotation symmetry breaking at the mean-field level, 
it appears that TZqq, V,PH, To t<->J, and T o TTZ are 
always good symmetries. After the breaking of the S z 
rotation the above symmetries might be modified; after 
the operation of the discrete symmetries it might require 
an additional S z rotation to restore the ground state. 
For example TZ§q and V are modified to e l7vSz o TZgo and 
e "rS* -p ^ ue orc i er . p or example, the angle of 

the S z rotation can be determined as follows. Assume 

e l6S - oTZ go \GS) = \GS) (C19) 

we obtain 

(GS\Sj\GS) = (GS\(e l9S * oTZ 9Q )Sl(e l9S * oTZgoy^GS) 
= e-^iGSlS^lGS) = -e- ie (GS\Sl\GS), (C20) 

hence 6 — ir. After some simple calculation it can be 
shown that in the antifcrromagnetic state PH and ToTTZ 
alone remain good symmetries. Finally, since To ]^{ 
transforms XY order parameter M to — M* (note how 
the nature of the AF order enters), the associated S z 
rotation must rotate — M* back to M. 

In summary, in the AF phase we find while T and TTZ 
are broken, e i7r5 * o^,e irf ' oV, ToTTZ and PH are good 
symmetries. This is consistent with a regular easy-plane 



AF ordered state with PH symmetry. Thus the AF phase 
portrayed by the slave-rotor+ gauge fluctuation theory is 
the regular AF phase. 
APPENDIX D: A MEAN-FIELD THEORY FOR 
THE AF PHASE AND THE ANALYSIS OF THE 
CORE STRUCTURE OF THE AF VORTICES 

In the main text we have started with a mean-field 
superconductor and show its vortices carries spin. The 
condensation of these vortices destroy the superconduc- 
tor and drive the system into an insulating easy-plane 
antiferromagnet. The vortices of this antiferromagnet 
carries charge and their condensation destroys the AF 
order and drives the system back to the superconducting 
phase. In section [C] we have presented a gauged slave 
rotor theory for this antiferromagnet. In this appendix 
we present a mean-field Hamiltonian for this phase and 
study its vortices. 

From Sec[C] we have shown the symmetries of the AF 
phase are charge-Z7(l), e i7rS * oTZ go , e i7rS * oV, PH, ToTTZ. 
These symmetries severely limit the possible quadratic 
mean- field Hamiltonians. The most general quadratic 
Hamiltonian can be written as combination of hopping 
terms (however in general the hopping does not have 
to conserve the spin). Demanding the above symme- 
tries it is possible to write down all symmetry allowed 
hopping terms. For example consider the nearest neigh- 
bor hopping. We shall show in the following that the 
only symmetry allowed such hopping term is of the form 
Xi(c\ a Cj C + h.c.) where xi us a rea l number. Let us 
first assume the up spin hopping amplitude is imaginary: 
iXc] 

t c iT + ^- c " Under PH transformation it becomes 

—iXciic^+h.c. = iACj^Cii+h.c. = — iAc^Cji+h.c. Hence 
the spin up and spin down hopping matrix elements must 
differ by a sign, i.e., ia(c\ (7 Cj a -—c^ a Ci lT )+h.c.. Because the 
sign difference between i — > j and j —> i, this hopping has 
a direction. Let us represent the direction in which the 
hopping amplitude is iX by an arrow. T o TTZ requires 
the arrows to be translation invariant, but this breaks 
the e mS * oR 90 . So the imaginary spin preserving nearest 
neighbor hopping is forbidden. Next let us consider the 
spin flipping nearest neighbor hopping M*c\^Cj[ + h.c. 

PH transforms it into M*CiicL+h.c. = — M*cLcjj,+/i.c. 
Thus the hopping also has a direction, i.e., the hopping 
amplitude for c\^Cji is the negative of that of c^Ci[. The 
T o TTZ requires the hopping to be translation invariant, 
but this again must break the e mS * o TZ%$ . So the spin 
flipping nearest neighbor hopping is also forbidden. This 
leaves the regular real-amplitude spin-preserving hopping 
as the only possibility. 

One can perform similar analysis for all possible 
quadratic terms: on-site terms, nearest neighbor terms, 
next nearest neighbor terms, etc. The following Hamil- 
tonian has the most general form up to second neighbor: 
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H AF = ^2Vi( M l4] C n + h - a ) + Xl ^2(4a C j" + h - C -) + X! tikivCivCkc + h - c -) + X] ^M^^l^i + c {] C *l) + h - C - 
i (ij) «ifc» «ifc» 

(Dl) 



where Uj has the staggered pattern as shown in FigQJa). 

Interestingly it is not possible to choose parameters so 
that Eq. ljDljl will have conic (i.e., Dirac-like) intersecting 
bands. If one set Mi = M 2 = 0, the two bands intersect 



quadratically at X point. Small M\ and/or M 2 open up 
a gap. In this limit the momentum space Hamiltonian 
around X is 
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where (fe 1? k 2 ) = k - X, M = M\ - 2M 2 and 

*l = ( C AT/C C STfe' C A|.fc' c Blk)- (D3) 

After a simple unitary rotation: 

- (; (D4) 

The form of the Hamiltonian becomes the same form 
as Eq. (|B9p in appendix B. Based on the index theorem 
obtained there, we conclude one vortex (antivortex) in 
M also has two mid-gap energy levels. If both levels 
are filled the vortex (anti-vortex) carries electric charge 
1 , and if both are empty the vortex (anti- vortex) carries 
charge —1. This is because the vortex with both levels 
filled and the vortex with both levels empty are related 
by PH and thus their charges differ by sign change. On 
the other hand by definition these two vortices' charges 
differ by 2. We thus prove the existence of the charge-±l 
vortex (anti-vortex) in the AF phase. 

APPENDIX E: THE MATRIX ELEMENT OF 

THE ORDER PARAMETER OPERATOR 
(VORTEX TUNNELING OPERATOR) IN SC 
AND AF PHASES 

The SC<->AF transition is realized by condensing vor- 
tices with non-trivial quantum number. The motivation 
of studying the vortex tunneling operator is to compute 
the order pattern in the vortex condensed phase. The 
calculation presented here provides a rigorous and sys- 
tematical way to determine the order pattern in the vor- 
tex condensed phase. As a first application, let us study 
the vortices in the SC phase to show that the vortex con- 
densed phase is AF. 



What is the order parameter in the SC vortex K^/f)! = 
7^ Ocondensed phase? Apparently both (Wf) and 
(^1) are non-zero and one may think that there are two 
spin-1/2 order parameters. However both \&-|- and ^ 
are non-local operators (since they annihilate vortices of 
the SC phase) hence can not serve as order parameter. 
However the combination ^j^x (which does not change 
vorticity but flips spin) is local, and acquires none zero 
value in the XY ordered state, hence can serve as an order 
parameter. To determine the nature (i.e. FM versus AF) 
of the XY order, we compute the quantum number of 
the ^j^j. operator. For example, if ^j^x changes sign 
under a 90 degree rotation (1Zgo) around the center of a 
plaquette , it is consistent with AF order. 

In order to determine the quantum number of the oper- 
ator ^j^x we consider the two vortex states |AVf) and 
|BV|) with the location of the vortices displaced from 
one another by one lattice spacing. However we tune 
the parameters so that the core size £ of the vortices is 
significantly larger than their separation. As shown in 
SecjB] there are two spin-up in-gap levels in the vortex 
core, and AV| and BVx differ by whether the two lev- 
els are filled or not. If 7^ and 7^ create Bogoliubov 
quasiparticle in these levels, then in the limit £ >> 1 
|AV T )=7j T 7^ T |BVx>. On the other hand from Eq. ([Aljl 

we know that the operator ^j^x sends BVx to AV|, i.e., 
it is the vortex tunneling operator: 

|AV t > = tff^BVj). (El) 

Therefore the quantum number of ^j^x ^ s simply the 
quantum number difference of |AVf) and |BVx). For ex- 
ample let us consider the rotation IZgo- If the vortex 
configuration is symmetric under 7^.90 we expect |AVf) 
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and |BV|) to be eigenstates of TZgo: (a) 

ft 90 |AV T ) =e 4e T|AV T ) 
n 90 \BV l ) = e i6 i\BV l ) (E2) 

To determine the quantum number difference, e 4 ^ - ^, 
we can calculate the following matrix element ratio 



(b) 



(AV T |S+|BV X ) 
(AV T |S+|BV X )' 



(E3) 



To see that let us consider j — lZgo(i) and thus Sj = 
90. In this case Eq. (|E3|) becomes 



"7^-go S~^7Z. 



(AV T |S+|BV 



(AVtI^ 



S+TZ 9Q IBV 



(E4) 



From the above it is clear that one can replace the S[ 
in Eq. (|E4[) by any spin-1 operator Oi and get the same 
result so long as (AV T |d t |BV|) ^ 0. 

In Eq. (|E4|) we have assumed that the vortex config- 
uration is 7?-9o-symmetric so that |AV|) and |BVj,) are 
eigenstates of Hqq. The easiest way to implement the 
90-degree rotation symmetry is to put in a vortex under 
open boundary condition. However this brings in edge 
states in the gap since the superconductor in question is a 
topological one. To avoid the complication of edge states 
it is better to use periodic boundary condition. However 
in that case one has to put in a vortex and an antivortex 
hence necessarily breaks the 90 degree rotation symme- 
try. Fortunately, it turns out that the ratio in Eq. (|E3|) is 
almost completely determined by the the wavefunctions 
of the mid-gap levels. The latter is localized in the core 
region and can only sense a region D around the center 
of the vortex. Thus if we separate the vortex and the 
antivortex by sufficient distance and make sure that the 
vortex configuration within region D is rotation symmet- 
ric, we expect the matrix element ratio would converge 
to the desired result e l ( e i~ e ^ in the following limit: (1) 
large vortex core size (2) large system size (D size). In 
practice we input the vortices by Jacobi theta function as 
discussed in detail in SeclEland the vortex configuration 
around the vortex center tends to rotational symmetric 
in the afore mentioned limit. 

In Fig[in]we present the result for (AV T 1 5+ 1 B V| ) as a 
function of i for a fixed vortex configuration. We find in 
the vortex core the matrix elements have staggered signs, 
which means the XY order is AF. 

We showed vE r j\f , x is an AF order parameter. Here we 
present some details of the numerics that we performed. 
In Fig[TT] we show the positions of the vortex and the 
anti- vortex centers on the torus. We use complex number 
Wi , W2 to represent the position of the vortex and anti- 
vortex. Given these positions, and using complex number 
z = x + iy to represent the position of the pairing bond 
center, the pairing amplitude is given by: 



A(z) = |A |e^ w (l - e- |z - tul|/ «)(l - e^-™ 2 ^), 



(E5) 
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FIG. 10: (AVf|S+|BVj.) as a function of i for sites close to 
the center of the vortex (labeled by ffi). We input a vortex 
and an anti- vortex separated by half system size on a 24x24x2 
('2' means 2 sites per unit cell) system on torus (details of the 
vortex configurations are given in SeclE}. Re(x), t and Im(x) 
are the parameters in Eq.([T}. To make the core size £ big, 
we input the pairing order parameter in the vicinity of the 
vortex core as Im(x)(l — e~ r ^)e' ,e , where r,9 are the polar 
coordinates with respect to the the center of the vortex. We 
used two sets of parameters: (a) large pairing and small t 
and (b) small pairing and large t. In limit (a) the vortex 
core is relatively small, while (b) the vortex core is large (for 
details see Sec|B]and|Ej. In both limit we see staggered signs, 
implying AF order. 




FIG. 11: (a) The real space unit cell of the square lattice 
model Eq.([TJ and the basis vectors a\,a,2. (b) In our numer- 
ical study, we put one vortex © and one anti-vortex sepa- 
rated by half system size on a sample with periodic boundary 
conditions. 



where £ <§; L is the size of the vortex core, and the phase 
4>{z) is given by: 



V2L 



27rlm(z 



L 2 



(E6) 



= arg (oi C^ 1 ' . < ■ ~ )^ r"-;-' . < 



And the Jacobi theta function is defined as: 

6i(z,q)= J2 (-ir^M"^ e {2n+1 ^ z . (E7) 

n— — 00 

One can show that e 1 ^^ is a periodic function on the 
LxL sample with a vortex and an anti- vortex at w\,W2- 
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We mention that £ determines the input vortex size, 
but it does not determine the size of the zero mode wave- 
function (which actually weakly dependent on £). Wc 
find that size of the zero mode wavefunction is mainly 
determined by the pairing |Ao| by numerical study. This 
behavior at least can be understood in the Dirac limit 
where the analytical solution of the zero mode wave- 
function is knowr*2£ to be ~ e~l Ao ' r . To provide a 
quantitative understanding of the size of the zero mode 
wavefunction ijj(r) we numerically fit the density 2 

(P _F 0) 2 

by Ae & . We find for the vortex configuration in 
Fig fTOT a) the zero mode size I = 1.5, and for Fig llO( b) 
I = 6.0. This finishes the discussion of the SC vortex 
tunneling amplitude. 

From now on we study the AF vortex tunneling ampli- 
tude. We adopt the mean- field Hamilton Eq. (|Dl[) from 
SeclDland only consider the on-site magnetization: 

Haf = J2^ M * c h Ci i + h - c -) + xJ2( c l c i° + h - c -) 



^ *»fc(o'ct r Cfc <T + h.c.) 
«ifc» 



(E8) 



From Sec|D]we know that there are two zero modes in the 
vortex core 771, 772 in the long wavelength limit. Therefore 
there are two charged vortex states |Vq=i) and |Vq=_i) 

are related by |Vq=i) = ^/i^^ I^Q=— 1) ■ O n the other hand 
because the z-bosons in Eq. (fT5]) is the relativistic bo- 
son fields of these vortices and anti- vortices, we conclude 
that: 



\Vq=x) =4z2\Vq=- 



(E9) 



This is because z[z2 does not change vorticity but create 
a pair of electrons. z\z-i is nothing but a cooper pair 
creation operator. In the (zi) 0, (22) 7^ phase z\zi 
acquires non-zero expectation value and it is the SC order 
parameter. What is the SC order pattern? 

To answer this question we perform the following ma- 
trix element ratio computation. 



{V Q=1 \{c\ ]C \ +{xy}l 



t „t 



(E10) 



where site i and site i + {x,y} form a n.n. bond. Simi- 
lar to the discussion in Sec MI) one can convince oneself 
that if the bond in the numerator and the bond in the 
denominator are related by 7?-9o, the 90 degree rotation 
around a plaquette center, then this ratio computes the 
TZqo quantum number of z\z-i operator. In Fig ll2l we see 
the staggered sign of the bond pairing matrix elements, 
which indicates the d-wave pairing symmetry in the AF 
vortex condensed phase. 
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FIG. 12: The AF vortex matrix elements 

(VQ=i\{c\^c\ +{xy}l - 4 i c i t +{a: ^ }T )|VQ = -i}. We input 
the vortex and anti-vortex in the same fashion as we did 
for SC vortex, and choose two sets of parameters: (a) large 
magnetization with small vortex size, and (b) small magne- 
tization with large vortex size. The staggered sign indicates 
that the pairing is d-wave in the AF vortex condensed phase. 



APPENDIX F: THE QUANTUM NUMBER OF 
THE MONOPOLE OPERATOR AND THE 
MAGNETIC ORDER PATTERN IN THE 
SLAVE-ROTOR THEORY 

In this section we determine the magnetic order pat- 
tern in the slave rotor theory. We first show the quantum 
number of the gauge monopole determines the magnetic 
order pattern, then we find a way to compute this quan- 
tum number. We find the magnetic order pattern is anti- 
ferromagnetic consistent with the conclusion in the main 
text. 

In the magnetic ordered state, the z a 's are absent at 
long distance/time. Dropping them from the slave-rotor 
gauge theory we arrive at a spinon theory interacting 
with a U(l) gauge field: 



iiatt = J2 fi^t + ib 0li - KAiyVi + [ J2 



i(6y-AfX)v..f+ 



Xijflfjv+ rr/ < 



r 



e i(6«-A,s a 3 ) f t 



«*.J»," 



fUia + h.C 



(Fl) 



Physically Eq. (jFlj) describes a half-filled (i.e., one parti- 
cle per site) system of fermions, or more precisely, a spin 
state. 

We now show that this theory describes an XY ordered 



phase. As discussed in the paragraph containing Eq. (|16p 
the monopole operator serves as the order parameter of 
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the magnetic order, i.e., we expect 



V, 



D 



(F2) 



where D is the spatial extent of the monopole (i.e., the 
region spanned by the inserted gauge flux). We expect 
to be the antiferromagnetic order parameter hence 



Vi 



(F3) 



where % is the staggered sign. In general there is an 
arbitrary global phase factor relating to J^ieD Vi^t > 
however this phase does not affect our determination of 
the the magnetic order. 

It is convenient to choose D and the inserted flux dis- 
tribution so that they do not break lattice symmetry of 
Eq. (|F1[) . The transformation properties of under 
the symmetry operation is what we referred to as the 
monopole quantum number. As discussed in appendix 
Othe symmetry of Eq. (IF1[) include rotation TZqq and/or 
ToTlZ. \ie %9i — IV? (ferromagnetic) Vp would be invari- 
ant under TZgct- On the other hand if e l6i = r\i (antifer- 
romagnetic) Vq would change sign under rotation. Thus 
the transformation property (or the quantum number) of 
crucially depends on the phase e ldi (or more precisely 
on the relative phase e l ( Si ~ s ^ for i 6 A or B sublattices. 

How to compute the relative phases? In the following 
we will show that 



f {GS, l-fluxIPS+PlCS', 0-flux) \ 
V (GS, l-flux|PS'+P|GS', 0-flux) / ' 



where | GS, 0-flux) is the spinon mean-field wavefunc- 
tion of Eq. @ in zero gauge flux and \GS, 1-flux) is the 
mean-field wavefunction in the background of a uniform 
flux integrated to one flux quantum. P is the opera- 
tor that executes the projection onto the Hilbert space 
of one-fermion-per-site. This formula of computing the 
monopole quantum numbers is firstly proposed by one 
of the authors in Ref. 38 to address the similar issue of 
the monopole quantum number in a U(l) Dirac spin liq- 
uid on a Kagome lattice. One way to understand this 
formula, is to consider the wavefunction of a monopole 
condensed phase. Clearly, it is a superposition of differ- 
ent fluxes \M) = P|0-flux) + aP| 1-flux) + a*P\-l flux). 
Hence the expectation value of the spin operator is given 
by (M\S+\M) ~ (0-flux| PS+P 1 1-flux). 

We have performed numerical calculation of the rela- 
tive phases in Eq. (1F4|) . Specifically, we fix i-site to be an 
arbitrary site on the square lattice and compute 6j — 8i as 
j varies through the sites of a finite lattice. (We employ 
the periodic boundary condition in this calculation.) The 
result for different lattice sizes are illustrated in Fig- P^)) . 
We find eWi-W is 1 if % and j are on the same sublattice 
and — 1 if i and j are on different sublattices. Moreover 
the result is independent of the value of x an d t. This 
proves that the magnetic order described by Eq. (|F1[) is 
anti-ferromagnetic. 
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FIG. 13: The numerical result for Eq. ||F4) for 4 x 4, 8 x 8, 12 x 
12 and 16 x 16 lattices under periodic boundary condition. 
We fix i to be the bottom left site and show the result as j 
goes through six different sites as shown. In the limit of large 
system size the magnitude of the ratio approaches one and 
the phases factor is 1 for i,j on the same sublattice and —1 for 
the opposite sublattices. In constructing this fiture we choose 
parameters $ = 0.67T, \t/\\ = 0.2 in Eq. ((2]). The relative 
phase factor are found to be independent of the parameter 
values. 



In the following we present a discussion of why Eq. (|F4j) 
is the right quantity to calculate. As discussed in ap- 
pendix [Clsvmmetrv operation is manifested in gauge the- 
ory modulo gauge transformation. In other words, any 
change of the /-Hamiltonian in Eq. (|Flj) induced by a 
transformation that can be undone by a gauge transfor- 
mation will preserve the gauge theory defined in Eq. (|F1[) . 
The collection of such compounded transformation form 
the so-call projective symmetry group (PSG) 3 —. How- 
ever, unlike the physical symmetry operations, the ele- 
ments of the PSG depends on which gauge the mean- 
field /-Hamiltonian is written. Certain gauge are par- 
ticularly convenient (in the sense that the accompanied 
gauge transformation is trivial) for some physical sym- 
metry operations (but not for others). However, so long 
as two gauges are related by a gauge transformation 
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where e^' = 1 



(F5) 



the corresponding PSG operation will yield the same re- 
sult when acting upon a spin state (i.e., a spinon state 
with one particle per site) . This degrees of freedom works 
in our favor when we try to determine the symmetry 
properties of spin state. For example, it allows us to 
choose the "best gauge" for each different symmetry op- 
eration. Thus in order to determine the symmetry prop- 
erty of a spin state (written in spinon variables) one can 
adapt the "best gauge" for each symmetry operations. 

Now we are ready to define the monopole quantum 
number. Let |W, 0-flux) and |\&, 1-flux) be the ground 
spin state (written in terms of the spinon variables) with 
zero and one quantum of background flux, respectively. 
(By definition |#, 1-flux) = Vt|#, 0-flux).) The quantum 
number of is determined by the relative transforma- 
tion properties of I*?, 1-flux) and |^, 0-flux) under the 
elements of PSG (hopefully under the best gauge). Let 
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FIG. 14: The in-gap landau levels when a single quantum 
of magnetic flux is added to Eq. ([2}. K and K* mark the 
momentum space location of the Dirac points when the second 
neighbor hopping is switched off. For up-spinons the Landau 
level are close to the top of the valence band, and for down 
spinons they are close to the bottom of the conduction band. 
In the depicted filling of Landau levels, the total S z = 1. 

us again use rotation IZgo as an example. If |'F, 0-flux) is 
invariant but \^f, 1-flux) acquires a phase factor e 1 ^ under 
the appropriate PSG operation 

ft 90 0-flux) = | tf, 0-flux) 

K 9Q \y, 1-flux) = e^|*, 1-flux), (F6) 

then by definition V' has a rotational quantum num- 
ber e^. Because PIGS', 0-flux) and P\GS, 1-flux) also 
describes state with one spinon per site, they also en- 
joy the best gauge freedom discussed earlier. It is also 
important to point out that although PIGS', 0-flux) and 
P | GS, 1-flux) arc not the exact states I 1 ?, 0-flux) and 
1^, 1-flux). So long as they have the same PSG trans- 
formation property as them, Eq. (|F4|) (in the thermody- 
namic limit) will yield the exact monopole quantum num- 
ber. Again using TZqq as an example we should have 

(GS, l-flux|PS+P|GS, 0-flux) 
(GS, l-flux|PS+ 9o(i) P|GS, 0-flux) 

(GS, l-flux|PS+P|GS, 0-flux) 
~ (GS, l-B.vx\PK^SfK9oP\GS, Mux) ' 

The advantage of P|GS, 0-flux) and P|GS, 1-flux) are 
that their PSG transformation properties can be be read 
off from those of \GS, 0-flux) and \GS, 1-flux). 

Even so determine the transformation properties of 
| GS, 0-flux) and | GS, 1-flux) can be quite involved. For 
example if we want to obtain the transformation prop- 
erty of | GS, 0-flux) under TZgo, one can choose a gauge 
(the best gauge) in which the the spinon hopping matrix 
elements are fully translation invariant. In this gauge we 
expect | GS, 0-flux) to be rotation symmetric. However 
to figure out the symmetry properties of |GS, 1-flux), we 
face a more challenging situation. For example it is im- 
possible to choose a gauge which preserves the rotation 
symmetry on torus (because the fluxes through the torus 
holes are gauge invariant and breaks rotation). Under 
open boundary condition one can choose the symmetric 
gauge A = \Bz x r which is rotation invariant. Under 



this (best) gauge it is possible to analytically study the 
transformation law of |GS, 1-flux) under rotation. More 
specifically, due to the topological nature of the spinon 
mean-field Hamiltonian, one can show that in the pres- 
ence of one quantum of uniform flux, there are two in-gap 
Landau levels near the top of the valence band for the 
up spinons and another two near the bottom of the con- 
duction band for the down spinons. This is illustrated in 
FiglHl Because these in-gap levels have definite angular 
momentum in the symmetric gauge an analytical study 
of the rotation quantum number is possible. In numerical 
study of Eq. (|F4[) one does not need to worry about any 
of the above; one simply computed in the ratio in the any 
gauge. We have done just that for the results presented 
in FigCTJ 

APPENDIX G: SOME SMALL SECTIONS 

1. Microscopic symmetries and action Eq. (|13[) 

As a double check we show in the following that the 
validity of action Eq. (|13[) is ensured by microscopic sym- 
metries. We should check that under symmetry transfor- 
mations (1) z\ — > Z2 (2) z — > z^ (which ensures theory to 
be relativistic), and (3) z\ and zi particles do not experi- 
ence any background magnetic field. In the AF phase the 
spin S z rotation is broken and charge-L T (l), PH, ToTIZ 
remains good symmetries. Due to the antiferromagnetic 
order, reflection V has to be followed by a 7r spin S z ro- 
tation, i.e., V — > e lnSz o V so that the ground state will 
be preserved. 

PH transforms a spin vortex into a vortex and flip the 
charge, therefore 

PH: Zl -> z 2 , z 2 -> z x (Gl) 

Under e l7rSz o V a vortex is transformed into an anti- 
vortex while preserving the charge, thus we have 

e™ s * oV: zi^> zl, z 2 -> z\ (G2) 

Eq. (|GllG2|) ensures condition (1) and (2). ToTTZ trans- 
forms a vortex into a vortex while preserving the charge: 

T o T1Z: zi — > zi Z2 — > Z2 (G3) 

Although z a transform trivially under T o T1Z, it is im- 
portant to bare in mind that because it is an anti-unitary 
transformation it changes sign of in Eq. (|13p . This en- 
sures condition (3). 

2. The edge modes of the topological 
superconductor 

Since the rotation Eq. ([3]) preserves the spectrum of the 
Hamiltonian, we can use the existence of gapless edge 
states in TBI to imply the same for the SC. Moreover 
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since the rotation preserves the spin, we conclude that, as 
in TBI, the edge modes in the superconductor must have 
opposite spin propagate in opposite direction. However 
unlike the edge mode in TBI, that of the superconductor 
do not carry definite charge. 

It is important to ask under what condition are these 
edge modes stable. In the TBI phase it is the charge 
and spin conservation law that guarantee the stability. 
Translated into the SC, spin conservation law remains 
the same but the charge conservation law becomes the 
conservation of J2i Vii^^l + h.c), does not seem to be 
a physical conservation law for a realistic system. Fortu- 
nately it turns out that the particle-hole symmetry plus 
the S z conservation are sufficient to guarantee the stabil- 
ity of the gapless edge modes in the SC phase. 

This is easily seen as follows. Particle hole symmetry 
is implemented by tp r — > i[i x tp r . For convenience, let us 
translate this into the spinon variables where it implies 
F r — > Z r i[i x Z\F r , which is F r — > iji z F r in a particular 
gauge. If the quadratic Hamiltonian of the F r is invariant 
under this transformation, then it should only depend on 
the identity and /i z . This implies it has a U(l) rotation 
symmetry corresponding to spinon number conservation. 
Combined with U(l) spin rotation invariance, this is suf- 
ficient to protect the edge modes. 



3. Physical Interpretation of the doublet of 
Chargeons 

Since the preceeding justification for introducing a 
charged boson doublet is rather formal, phrased in terms 
of gauge theories and symmetries, we provide here a more 
physical motivation for the appearance of the pair of 
fields z\, Z2- Consider an insulating state composed of 
local moments described by a fermionic mean field the- 
ory as in !25l Consider introducing charge fluctuations by 
including electron excitations c\ . In a single band model, 
these electrons or holes will form a singlet with the spin 
on that site. Hence, if we consider low energy charge ex- 
citations which are spin singlets, then one can construct 
two different operators b\ = c£/ CT and b\ = ^aa'c\j\, 1 
where e aa i is the antisymmetric symbol and a sum on 
spin indices is assumed. Keeping in mind the constraint 
|2"21 one can readily show that these bosonic operators sat- 
isfy: 

c T = 6i/t - /|4 
Crl = hf ri + f\b\ (G4) 

clearly these are very similar to the z\, Z2 fields. More- 
over, one can check that the on site Cooper pair operator: 
£aa'ct c lr' = Hence, one can consider these opera- 

tors as a 'Schwinger boson' representation of the Cooper 
pair. 



APPENDIX H: XY^SC DIRECT TRANSITION 
ON THE HONEYCOMB LATTICE 

From the vortex core zero mode study (see SecHK ITVl) 
we found there is a direct transition from AF phase to d- 
wave SC phase on square lattice. However that approach 
does not seem to have a predicting power. For example 
it is unclear how to generalize the square lattice results 
to a different lattice. On the other hand the slave-rotor 
theory (see Sec|V} can reproduce all those results in a 
systematical way, and it is obviously can be applied to 
other lattices. In this section we will apply the slave- 
rotor theory on honeycomb lattice, following exactly the 
same route as in SecfV] and discuss its XY<-^SC direct 
transition. 

We again firstly present the minimal requirement 
of symmetries for microscopic models containing the 
XY<->SC direct transition. In summary the symmetries 
are: charge-£/(l), S z — U(l), PH (particle- hole), V (re- 
flection of the mirror along the diagonal direction of 
the hexagonal plaquette) and TTZ (TTZ is time-reversal) . 
These transformations affect the electron operators as 
follows: 
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To describe the XY magnetic ordered phase on the 
honeycomb lattice by slave-rotor theory, we again assume 
that the bosonic chargon z are fully gapped, and the 
fermionic spinon / band structure to be TBI 22 : 

/// /; ' = -t E flfj* + * E A <,< (/:•//• - flifn) 

m did)) 

+ h.c, (H2) 

where \j = A tj x dl ? z = UijX, and dji,du are the two 

bonds connecting site j and i, Vij = ±1. This band struc- 
ture is the same as the Kane-Mele spin-hall insulator— if 
one replaces the spinon operator / by electron operator 
c. But the difference is that here the spinon / automat- 
ically couples with a dynamical 2+1D ?7(1) gauge field. 
Due to the "spin-hall" effect of the TBI band structure, 
we learn that the photon mode of the U(l) gauge field is 
actually the Goldstone mode of the XY magnet. Similar 
to the square lattice case, the low energy effective theory 
of the XY magnet has the form of Eq. (fl3|) . One can show 
that the microscopic symmetries Eq. (|H1[) ensure the cor- 
rectness of this low energy effective theory. 

The issue is to determine the XY and SC order pat- 
terns. We first study the XY order pattern. Is it ferro- 
magnetic or anti-ferromagnetic? To answer this question 
we again perform two calculations both showing the or- 
der is AF. First we computed the matrix element ratio 
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4 by 4 results 



8 by 8 results 



FIG. 15: We present the numerical result of the ratio of ma- 



trix elements 



(gs,1-&ox\ps+ p\gs,0-&ox) 



(GS ,l-fiux|PS+p|GS,(Miux> for the P ra J ecte d TBI 
wavefunction on 4 by 4 and 8 by 8 unit cell lattices on torus. 
We always fix site i to be the left-bottom one and show the 
result of site j on figure. The P\GS, 1-flux) breaks physical 
translations and rotations since there are fluxes in the torus 
holes, and thus the magnitude of this ratio is not exactly 
one. But in the thermodynamic limit where the boundary 
condition is irrelevant, the magnitude of the ratio should ap- 
proaches one, which seems to happen if one compares the 4 
by 4 result and 8 by 8 result. The relative phases of the 
matrix element is exactly for same sublattice and n for op- 
posite sublattices. We choose hopping parameters X/t — 0.1 
in these calculations, but the relative phases are found to be 
independent of choice of the parameters's value. 





(b) 



FIG. 16: (a) The honeycomb lattice in real space. The two 
sublattices are labeled by A and B. a\ and 0,2 are the two 
unit cell basis vectors, (b) The Brillouin zone of the honey- 
comb lattice. The vectors 61 and 62 are the two basis of the 



reciprocal lattice satisfying 



bj = Sij. 



to study the quantum number of monopoles of the U(l) 
gauge field (see SecJF]). We find the magnetic order is 
anti-ferromagnetic and in Fig. [15] we presented the ma- 
trix element ratio results. 

Secondly we directly measure the spin-spin correlation 
function of the projected TBI wavefunction. We label 
the two sites in one honeycomb unit cell by A and B as 
shown in FigOU and define |* ) = P\GS, 0-flux). The 
correlation functions that we measure are: 



C AA (r) = (Sl(r)S%(0)- 
C AB (r) = {S%(r)S%(0) 



S y A (r)S» A (0)U o , 

s y B (?)S y A (o)h . 



(H3) 



Because the projected wavefunction is the simplest way 
to mimic the effect of the gauge fluctuations (i.e., the 
one-fermion-per-site constraint), we expect that this cor- 
relation function to show AF order. 




FIG. 17: For 8 by 8 unit cells, we fix hopping parameters 
t = 1, A = 0.1 and plot Fourier transform of the real space 
spin spin correlation function {S x (r)S x {0) + S y (r)S y (0)} of 
the projected wavefunction in the hexagonal Brillouin zone 
of the honeycomb lattice, (a) is the correlation function of 
the same sublattice C AA . (b) is the correlation function of 
different sublattices C AB . The large q — positive peak for 
C AA and negative peak for C AB is a strong evidence that this 
projected wavefunction is AF ordered. 



We firstly present the projected wavefunction result for 
a fixed lattice size, namely 8 by 8 unit cells (see Fig HT)) . 
The momentum space peak at q = for Fourier trans- 
form of the XY spin correlation (positive for C AA and 
negative for C AB ) is a strong indication that the mag- 
netic order is AF. 

Next we present the finite size scaling study. When 
X/t = 0, the projected wavefunction has two Dirac nodes 
and is expected to show power law scalings of the spin 
correlation function. We choose periodic boundary con- 
dition for difference system sizes L = 8, 10, 14, 16, 20, 22 
(Note that 6N x 6N lattice has nodal fermions, namely 
fcrmion modes right at the nodes. Thus we avoid those 
lattice sizes.), and compute the spin correlation functions 
C AA and C AB at the half way of the system sizes. 

We checked that the projected wavefunctions that we 
used respect the full lattice space group symmetry. By 
the lattice space group symmetry, we find the three 
spin correlations at half system sizes for C AA , namely 
C AA {^), C AA {^) and C AA ( L{ai + a2} ), are related 
by symmetry transformation and thus identical. Here a\ 
and a 2 are the real space unit cell basis of the honeycomb 
lattice as shown in FigUB] On the other hand for C AB 
we find C AB (^) and C AB ( Mzi+^A ) are identical, but 
C AB {^ L ) are different from them. Therefore we present 
the results of C AA {±f-), C AB (^) and C AB (^f-) in 
FigHHa). (Note that C AA data are positive and C AB 
data are negative, and we present the absolute values 
only.) We find 



X/t = 0: C(r) cx 



f|1.89(2) 



(H4) 



Now we turn on A. Based on the low energy effective 
theory (see text) we know the ground state should be 
long-range XY spin ordered. But by studying the pro- 
jected wavefunction, we find the spin correlation func- 
tion still have power law scaling. The scaling exponent, 
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FIG. 18: We present the log-log plot of the finite size scaling 
of the spin-spin correlation functions for parameters t = 1 and 
(a) A = and (b) A = 0.05 and 0.1. From (a), we find when 
X/t = the scaling power law is \C\ oc x g 9(2) . From (b) 
we also find power-law scalings for X/t — 0.05 and X/t — 0.1, 
and the exponents for these two cases are the same up to 
statistical error, which is C oc t- oo(a) • The dashed line in (b) 
is the plot of C = 2^. 

however, is strongly modified from the A = value (see 
FigHHb)). We find 

X/t = 0.05,0.1: C(r) oc (H5) 

This means first the projected wavefunction fails to pro- 
duce the true long-range correlation, and second the AF 
correlation is strongly enhanced with A ^ projected 
wavefunction compared with the A = one. 

Now we study the SC order pattern. This is straight- 
forward because we should simply condense z\ and zi 
together. Deep in the SC phase the easy-plane conden- 
sation of z-boson can be achieved by doing replacement 
z x = -^e iSl and z 2 = -^e lf>2 in Eq. (|2"5]l: 

c it = ^(/ite <fll -»7i/le- <fla ) (H6) 
C i = i(/iie tfl +iji/le-*). (H7) 



We can express / fermions in terms of the electron oper- 
ator c: 

/it^tete-^+^e-*") (H8) 
f il = 1=^-^-7,4^) (H9) 

Plugging these into Eq. (|H2p . the superconductor Hamil- 
tonian is: 

H sc =t c\ a c ja +i\Y, ^me l(<>2 ~ ei) {c^c 3l + c tl c n ) 
(id) ((id)) 

+ h.c. (H10) 

If we absorb e l ^ 2_ei ' into A so that A is complex, the 
full n.n.n. pairings are 

(-fuij [iA(c iT c^ + c^cjf) + iX*{c\4 3 i + c \i c ]])^] 

(Hll) 

these are /-wave triplet pairings with S z = 0. Note 
that this /-wave SC is fully gapped (its spectrum is the 
same as TBI). Therefore we show that there is a direct 
AF<->/-wave SC transition on honeycomb lattice, whose 
low energy effective theory is again the easy-plane limit 
ofNCCP 1 . 

Finally we remark that there are zero modes in the 
vortex cores of the AF and f-wave SC phases by solving 
the BdG equations. We already presented the mean-field 
Hamiltonian of the SC phase in Eq. (|H10[ ). To write down 
the AF mean-field Hamiltonian one should study which 
terms are allowed by symmetry (in AF phase, the good 
symmetries are charge-J7(l), PH, V and e l7rSz o Til). 
The simplest AF mean-field Hamiltonian consistent with 
symmetries is: 

H XY = tJ24 a c ja + ffc^(M*cJ TCa + Mc\ [Cl] ). 
(ij) « 

(H12) 

where M is the on-site magnetization. 

When the order parameters are small both Eq. (|H10[ ) 
and Eq.(|Hl2j are in the Dirac limit. The AF side BdG 
equation has already been solved by Herbut^i. And one 
can show that the SC side BdG equation turns out to 
be the same as the AF side BdG equation after appro- 
priate rotations. One can show that the whole analysis 
on the square lattice, namely zero modes^vortex carry 
quantum numbers-direct transition, also carries through 
on the honeycomb lattice and reproduce the same low 
energy effective theory as the one given by slave-rotor 
theory. 
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